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SUMMARY 

Meteorological  icing  data  obtained  in  flight  in  the  United  States 
are  analyzed  statistically  and  methods  are  developed  for  the  determina¬ 
tion  of:  (l)  the  various  simultaneous  combinations  of  the  three  basic 
icing  parameters  (liquid -water  content,  drop  diameter,  and  temperature) 
which  would  have  equal  probability  of  being  exceeded  in  flight  in  any 
random  icing  encounter;  and  (2)  the  probability  of  exceeding  any  speci¬ 
fied  group  of  values  of  liquid-water  content  associated  simultaneously 
with  temperature  and  drop-diameter  values  lying  within  specified  ranges. 
The  methods  are  particularly  useful  in  the  design  of  anti-icing  equip¬ 
ment  intended  to  operate  through  the  United  States,  to  define  simulta¬ 
neous  combinations  of  the  meteorological  variables  which  could  be 
encountered,  and  to  ascertain  the  effectiveness  of  the  equipment  in 
withstanding  the  natural  icing  conditions  to  which  it  may  be  subjected. 
In  addition,  a  mathematical  basis  is  provided  for  the  future  statistical 
analysis  of  meteorological  icing  data  that  might  be  obtained  throughout 
the  world. 


INTRODUCTION 


The  program  of  research  in  aircraft  ice  prevention  which  has  been 
conducted  by  the  NACA  during  the  past  several  years  has  been  directed 
primarily  toward  the  development  of  practical  methods  for  the  design  of 
thermal  ice -prevention  equipment  for  various  airplane  components.  Since 
a  rational  ice-prevention  design  requires  a  knowledge  of  the  physical 
characteristics  of  icing  conditions,  an  important  phase  of  the  research 
program  has  been  an  investigation  of  the  meteorological  conditions  con¬ 
ducive  to  icing. 

The  severity  of  an  encounter  with  icing. conditions  is  determined 
principally  by  four  factors;  namely,  the  liquid-water  content,  the 
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diameter  of  the  drops,  the  temperature,  and  the  horizontal  extent  of  the 
conditions.  The  meteorological  investigation  has,  therefore,  been  con¬ 
cerned  with  obtaining  measurements  of  these  quantities  in  a  wide  variety  w 

of  natural  icing  conditions  in  order  to  establish  the  range  and  rela¬ 
tive  frequency  of  occurrence  of  various  values,  and  combinations  of 
values,  of  these  quantities. 

During  the  period  1945  through  1948,  a  considerable  amount  of  data 
was  collected  by  the  Ames  and  Lewis  Laboratories.  These  results,  and 
discussions  of  various  phases  of  the  investigation,  have  been  reported 
in  references  1  through  4.  All  of  these  results,  together  with  similar 
data  from  observations  made  by  other  organizations,  were  used  as  a  basis 
for  a  listing  of  estimated  maximum  icing  conditions  recommended  for  con¬ 
sideration  in  the  design  of  anti-icing  equipment  (reference  5).  The 
probable  maximum  values  listed  in  reference  5  vere  estimated  on  the 
basis  of  what  was  considered  a  reasonable  extrapolation  from  a  limited 
amount  of  data,  and  the  concept  of  the  probability  of  encountering 
various  simultaneous  combinations  of  values  entered  only  in  a  subjective 
and  qualitative  way. 

In  this  report,  the  available  data  obtained  in  icing  flights  in  the 
United  States  are  analyzed  statistically  and  the  results  are  presented  ^ 

graphically  in  two  forms.  In  the  first  form,  the  results  are  plotted  to 
show  the  various  combinations  of  liquid-water  content,  drop  diameter, 
and  ambient-air  temperature  which  have  equal  probability  of  being 
exceeded.  In  the  second  form,  plots  are  presented  to  determine  the 
probability  of  exceeding  any  specified  value  of  liquid-water  content 
under  the  condition  that  the  value  is  associated  simultaneously  with 
values  of  temperature  and  drop  diameter  lying  within  specified  intervals. 

Both  plots  are  based  on  preselected  values  of  horizontal  extent;  however, 
methods  for  estimating  values  applicable  to  other  horizontal  extents  are 
also  presented. 

Although  the  results  of  this  report  are  directly  applicable  only  to 
the  United  States,  an  estimation  of  the  conditions  prevailing  in  other 
localities  of  similar  climatic  conditions  can  be  made  based  upon  the 
results  presented  herein.  As  meteorological  icing  data  are  accumulated 
for  regions  outside  the  United  States,  the  methods  of  this  report  should 
provide  a  framework  for  placing  the  data  on  a  statistical  basis  that  is 
not  limited  in  scope  to  the  United  States. 


SYMBOLS 


The  following  symbols  are  used  throughout  this  report: 

ft 

A,B,C  designation  applied  to  three  events  occurring  simultaneously 

but  not  necessarily  independently 
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c 

chord,  length,  feet 

D 

mean-effective  drop  diameter,  microns 

Dcr 

the  smallest  drop  diameter  for  which  drops  will  impinge  upon 
a  thermal  element 

AD 

interval  of  drop  diameter,  microns 

e 

Naperian  base  of  logarithms 

collection  efficiency  of  an  airfoil  section,  percent 

F 

factor  by  which  liquid-water-content  values  are  multiplied 
to  make  them  applicable  to  other  distances  of  cloud  hori¬ 
zontal  extent,  dimensionless 

% 

Mt 

weight  rate  of  water  intercepted,  pounds  per  hour  per  foot 
of  airfoil  span 

* 

n 

an  integer 

0 

an  icing  condition  considered  to  be  defined  by  specifications 
of  the  three  variables,  liquid-water  content,  temperature 
depression  below  freezing,  and  drop  diameter 

P 

probability  of  overloading  any  given  element  of  a  thermal 
anti-icing  system 

PA 

probability  of  the  occurrence  of  event  A 

PABC 

probability  of  the  simultaneous  occurrence  of  three 
events,  A,  B,  and  C 

pb(a) 

probability  of  the  occurrence  of  event  B,  under  the  condi¬ 
tion  that  event  A  will  occur 

PC(A,B) 

probability  of  the  occurrence  of  event  C  under  the  condi¬ 
tion  that  both  events  A  and  B  will  occur 

Pe 

probability  that  any  random  icing  encounter  will  be  char¬ 
acterized  by  a  combination  of  values  of  the  variables 
(liquid-water  content,  temperature  depression  below  freez¬ 
ing,  and  drop  diameter),  in  which  each  variable  must,  res- 
.  pectively,  equal  or  exceed  a  specified  value  of  that 
variable 

PG 

probability  defined  by  the  Gumbel  probability  distribution 
and  expressed  by  equation  (Al) 
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AP± 

probability  that,  for  a  random  icing  encounter,  the  following 
three  conditions  will  simultaneously  be  fulfilled: 

* 

(1)  liquid-water  content  exceeding  a  specified  value; 

(2)  temperature  within  a  specified  interval;  and  (3)  drop 
diameter  within  a  specified  interval 

W 

PAD 

probability  that  drop  diameter  will  have  a  value  lying  within 
a  specified  range  (AD) 

PAT 

probability  that  temperature  will  have  a  value  lying  within 
a  specified  range  (AT) 

PWi(AT) 

probability  that  liquid-water  content  will  be  equal  to  or 
greater  than  ¥j_  under  the  condition  that  the  temperature 
lies  within  a  specified  range  (AT) 

Q 

probability  defined  by  the  probability  function  Q(T' )  for 
a  fixed  value  of  T' 

Q(T') 

probability  function,  determined  by  the  variation  with  tem¬ 
perature  depression  below  freezing  for  the  most  severe 
icing  conditions,  of  the  percent  of  icing  encounters  in 

which  the  temperature  depression  below  freezing  is  less 
than  certain  specified  values 

V 

R 

probability  that  the  maximum  liquid-water  content  in  an  icing 
encounter  is  equal  to  or  greater  than  ¥,  under  the  condi¬ 
tion  that  the  temperature  depression  below  freezing  is 
equal  to  or  greater  than  T’ 

R(¥,T') 

probability  function  expressing  R  in  terms  of  W  and  T' 

R' 

probability  that  the  liquid-water  content  associated  with 
the  maximum  value  of  drop  diameter  in  an  icing  encounter  is 
equal  to  or  greater  than  W,  under  the  conditions  that  the 
temperature  depression  below  freezing  is  equal  to  or 
greater  than  T'  and  the  maximum  value  of  drop  diameter 
is  equal  to  or  greater  than  D 

R*(¥,D,T») 

probability  function  expressing  R'  in  terms  of  W,  D,  and  T' 

S 

probability  that  the  drop  diameter  is  equal  to  or  greater 
than  D,  under  the  conditions  that  the  temperature  depres¬ 
sion  below  freezing  is  equal  to  or  greater  than  T'  and 
the  liquid-water  content  is  equal  to  or  greater  than  ¥ 

S(D,W,T») 

probability  function  expressing  S  in  terms  of  D,  ¥,  and  T' 

S’ 

probability  that  the  maximum  value  of  drop  diameter  in  an 
icing  encounter  is  equal  to  or  greater  than  D,  under  the 

condition  that  the  temperature  depression  below  freezing 
is  equal  to  or  greater  than  T' 
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S'(D,T' ) 
S'  '(D) 

t 

T 

T' 

Tcr' 

AT 

AT' 

u 

V 

W 

Wo  .05 

Wo  .63 

AW 

X 


probability  function  expressing  S'  in  terms  of  D  and  T* 

composite  distribution  function  consisting  of  parts  of  each 
of  the  functions  S(D,0,0')  and  S'(D,0) 

airfoil  thickness,  percent  of  chord  length 

temperature  existing  in  an  icing  condition,  °F 

temperature  depression  below  freezing,  equal  to  32 -T,  F° 

temperature  depression  below  freezing  at  which  the  heat  out¬ 
put  of  the  thermal  element  can  bring  the  surface  tempera¬ 
ture  of  the  element  just  to  freezing  temperature ,  F° 

interval  of  temperature,  °F 

interval  of  temperature  depression  below  freezing,  F° 

a  constant  which  denotes  the  mode  in  the  Gumbel  probability 
equation 

airspeed,  miles  per  hour 

liquid-water  content,  grams  per  cubic  meter 

value  of  liquid-water  content  corresponding  to  a  probability 
value  of  0.05  defined  by  the  Gumbel  probability  equation 

value  of  liquid-water  content  corresponding  to  a  probability 
value  of  O.63  defined  by  the  Gumbel  probability  equation 

interval  of  liquid-water  content,  grams  per  cubic  meter 

independent  variable  in  the  Gumbel  probability  equation 


a  a  parameter  which  depicts  the  concentration  of  a  frequency 

distribution  about  the  mode  in  the  Gumbel  distribution 
equation 


Subscripts 


1 


any  selected  icing  condition 


i 


a  representative  value  of  a  variable 


ANALYSIS 


The  data  used  in  this  analysis  consist  of  1038  sets  of  values  of 
liquid-water  content,  drop  diameter,  and  temperature  obtained  during 
252  encounters  with  icing  conditions.  Values  of  liquid-water  content 
and  drop  diameter  were  measured  by  means  of  the  rotating  multi cylinder 
method.  All  data  available  at  the  time  of  preparation  of  this  report 
were  used  in  the  analysis,  including  (l)  data  obtained  by  the  Ames  and 
Lewis  Laboratories  and  published  in  references  1,  2,  3>  and  6; 

(2)  data  contained  in  the  monthly  reports  of  the  Air  Materiel  Command 
Aeronautical  Ice  Research  Laboratory,  the  major  portion  of  which  is  pre- 
sented  in  reference  7j  and  (3)  data  obtained  by  United  Air  Lines  (refer- 
ence  8)  and  American  Airlines  (reference  9)  during  icing  flights  with 
Douglas  DC -6  type  airplanes. 

For  the  purposes  of  this  analysis,  the  basic  unit  of  data  is  the 
icing  encounter  which  consists  of  all  measurements  made  during  flight 
through  a  single  area  of  continuous  or  intermittent  icing  conditions. 

This  area  was  separated  from  other  icing  conditions3  by  relatively 
large  areas  in  which  no  icing  was  observed.  The  problem  considered 
herein  is  to  determine  the  probability  that  an  icing  encounter  chosen 
at  random  will  include  conditions  of  a  given  severity  averaged  over  a 
certain  distance.  Since  most  of  the  rotating-cylinder  observations 
represent  averages  over  distances  of  about  3  miles  in  cumulus  clouds 
and  10  miles  in  layer  clouds,  these  distances  are  regarded  for  purposes 
of  data  reduction  as  standard  values  of  horizontal  extent,  applicable 
to  data  from  the  two  principal  cloud  types.  The  probabilities  obtained 
on  this  basis  may  be  applied  to  other  distances  by  using  the  data  con¬ 
cerning  the  relation  between  horizontal  extent  and  average  liquid-water 
content  obtained  from  continuous  records  of  the  rotating-disk  icing-rate 
meter.  (See  reference  3.)  It  has  been  assumed  that  the  greatest  values 
of  liquid-water  content  and  drop  diameter  measured  by  the  rotating  cyl¬ 
inders  during  an  icing  encounter  represent  maximum  values  averaged  over 
the  standard  distances  of  3  or  10  miles.  This  assumption  is  probably 
approximately  true,  since  a  particular  effort  was  made  to  obtain 
rotating-cylinder  data  during  the  periods  of  most  rapid  ice  formation. 

Since  the  variations  of  temperature  during  a  single  icing  encounter 
are  usually  not  very  great,  the  most  severe  icing  conditions  measured 
during  a  particular  icing  encounter  then  will  be  either  the  greatest 
value  of  liquid-water  content  and  the  corresponding  value  of  drop  diameter, 
or  possibly,  the  greatest  value  of  drop  diameter  and  the  corresponding 
value  of  liquid-water  content.  Both  of  these  possibilities  for  a  severe 


3The  term  "icing  condition"  as  used  herein  denotes  a  state  of  the  atmos¬ 
phere  defined  by  a  set  of  values  of  temperature,  liquid-water  content, 
drop  diameter,  and  pressure  altitude  in  which  the  temperature  is  below 
freezing  and  the  liquid-water  content  is  greater  than  zero. _ 
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icing  condition  are  considered  in  the  probability  analysis  which,  for 
the  sake  of  convenience,  is  developed  in  three  phases. 

The  first  phase  of  the  analysis  is  the  classification  of  the  avail¬ 
able  data  according  to  cloud  type  and  geographical  location;  the  second 
phase  is  the  determination  of  the  particular  combinations  of  liquid- 
water  content,  drop  diameter,  and  ambient-air  temperature,  for  a  given 
horizontal  cloud  extent,  that  have  the  same  probability  of  being  simul¬ 
taneously  equaled  or  exceeded  in  a  random  icing  encounter;  and  the  third 
phase  is  the  determination  of  the  probability  of  exceeding  any  specified 
group  of  values  of  liquid-water  content  associated  simultaneously  with 
temperature  and  drop-diameter  values  lying  within  specified  ranges. 


Classification  of  Data  According  to  Cloud  Type 
and  Geographical  Location 


For  classification  purposes  in  this  report,  clouds  will  be  divided 
into  two  classes,  cumulus  clouds  and  layer  clouds;  and  the  area  of  the 
United  States  will  be  divided  into  three  regions,  the  Pacific  coast 
region,  the  plateau  region,  and  eastern  United  StateB  with  boundaries 
as  indicated  in  figure  1.  This  system  of  classification  divides  the 
observations  into  six  cases  as  indicated  in  table  I,  which  includes 
information  relating  to  the  number  of  icing  encounters  and  the  number 
of  individual  measurements  of  liquid-water  content  and  drop  diameter  for 
each  group.  The  two  cases  for  which  the  greatest  amount  of  data  are 
available,  eastern  layer  clouds  (case  5)  and  Pacific  coast  cumulus 
clouds  (case  2),  are  treated  in  considerable  detail.  The  other  cases 
are  treated  in  as  much  detail  as  the  data  permit,  with  the  exception  of 
case  6  (cumulus  clouds  in  eastern  United  States)  which  is  omitted  because 
of  insufficient  data. 


Determination  of  Particular  Combinations  of  Liquid-Water 
Content,  Drop  Diameter,  and  Temperature  Depression 
Having  the  Same  Probability  of  Being  Simul¬ 
taneously  Equaled  or  Exceeded  in  a 
Random  Icing  Encounter 


Exceedance  probability.-  The  three  fundamental  icing  variables  used 
to  define  an  arbitrary  icing  condition  are  liquid-water  content  W, 
drop  diameter  D,  and  temperature  T.  For  purposes  of  probability  anal¬ 
ysis,  however,  it  is  more  convenient  to  refer  temperature  to  the  freezing 
point  and  to  use  temperature  depression  below  freezing  T’  as  a  variable 
to  substitute  for  the  temperature  T.  Thus,  any  particular  icing  condi¬ 
tion  can  be  represented  by  selected  values  of  the  variables  Wx,  D-l, 
and  Tir.  The  probability  that  any  random  icing  encounter  of  standard 
extent  will  provide  simultaneous  values  of  each  of  the  icing  variables 
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equal  to  or  greater  than  each  of  the  selected  values  (i.e. ,  W>  ,  D>Di; 
and  will  be  termed  the  "exceedance  probability"  Pe  correspond¬ 

ing  to  the  icing  condition  Wi,  Di,  and  Ti*. 

The  exceedance  probability  represents  the  probability  of  the  simul¬ 
taneous  occurrence  of  three  events  which  are  not  necessarily  independent. 
In  order  to  evaluate  Pe,  use  is  made  of  the  theorem  of  compound  proba¬ 
bility  (reference  10)  which  states  that  the  probability  of  the  simulta¬ 
neous  occurrence  of  three  events,  A,  B,  and  C,  is  equal  to  the  probability 
of  the  occurrence  of  A  times  the  conditioned  probability  of  the  occur¬ 
rence  of  B,  under  the  condition  that  A  will  occur,  times  the  condi¬ 
tioned  probability  of  the  occurrence  of  C,  under  the  conditions  that 
both  A  and  B  will  occur.  This  theorem  can  be  expressed  symbolically 
as  follows: 


i>ABC  =  PA  pb(a)  pc(a,b)  (1) 

If  it  is  considered  that  the  probability  of  exceeding  certain  values 
of  T',  W,  and  D  is  analogous  to  the  probability  of  the  occurrence  of 
events  A,  B,  and  C  expressed  in  equation  (l),  then  the  exceedance 
probability  Pe  may  be  represented  by  the  product  of  three  probability 
factors  Q,  R,  and  S,  according  to  the  equation 

Pe  =  QRS  (2) 

In  applying  equation  (2)  to  the  meteorological  variables,  the  prob¬ 
ability  Q  may  be  considered  as  a  function  of  T'  alone,  R  a  function 
of  W  and  T',  and  S  a  function  of  D,  T’,  and  W.  The  functional  rela¬ 
tions  between  Q,  R,  and  S,  and  D,  T',  and  W  are  denoted,  respectively, 
by  the  expressions  Q(T’),  R(W,T')  and  S(D,T',W).  The  meaning  of  each 
expression  is  as  follows: 

Q(T')  a  relation  that  defines  the  probability  Q  in  terms  of  T* 

and  expresses  the  probability  for  a  given  value  of  T’ 
of  a  random  icing  encounter  having  a  temperature  depres¬ 
sion  below  freezing  equal  to  or  greater  than  T* 

R(W,T’)  a  relation  that  defines  the  probability  R  in  terms  of  W 

when  T’  is  specified  and  expresses  the  probability  for 
a  given  value  of  W  of  a  random  icing  encounter  having 
a  liquid-water  content  equal  to  or  greater  than  W,  under 
the  condition  that  the  temperature  depression  below  freez¬ 
ing  is  equal  to  or  greater  than  T* 

S(D,W,T* )  a  relation  that  defines  the  probability  S,  in  terms  of  D 
when  W  and  T'  are  specified  and  expresses  the  probabil¬ 
ity  for  a  given  value  of  D  of  a  random  icing  encounter 
having  a  value  of  drop  diameter  equal  to  or  greater 
than  D,  under  the  condition  that  the  liquid -water  content 
is  equal  to  or  greater  than  W  and  the  temperature 
depression  below  freezing  is  equal  to  or  greater  than  T* 
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The  physical  significance  of  the  term  exceedance  probability  is 
best  illustrated  by  graphical  means.  Figure  2  presents  a  three- 
dimensional  plot  with  the  variables  ¥ ,  D,  and  T*  as  the  axes.  Any 
icing  condition  of  standard  duration  defined  by  selected  values  of  these 
three  variables,  such  as  Dx,  and  Tx ' ,  appears  as  a  single  point  CK 

on  the  plot.  Corresponding  to  the  point  0^,  there  is  an  exceedance 
probability  P6l  which  expresses  the  probability  that  a  random  icing 
encounter  will  provide  simultaneous  values  of  each  of  the  icing  vari¬ 
ables  equal  to  or  greater  than  each  of  the  selected  values,  ¥x,  D,, 
and  ly.  According  to  equation  (2),  the  value  of  Pe  would  be  given 
by  the  equation  -1 


Pe1  — 

Thus,  when  Pe  is  regarded  as  a  constant,  Pe  ,  equation  (2)  defines  a 
surface  passing  through  02,  specified  by  all  the  combinations  of  ¥,  T', 
and  D  which  have  the  same  probability  of  being  simultaneously  equaled 
or  exceeded.  This  surface  is  shown  as  the  equiprobability  surface  in 
figure  2,  and  represents  the  locus  of  all  icing  conditions  which  have 
the  same  probability  of  being  equaled  or  exceeded.  The  quantitative 
value  of  exceedance  probability  applying  to  individual  points  lying  on 
the  surface  is  expressed  by  equation  (2). 

It  should  be  noted  that  since  Q,  R,  and  S  represent  probabilities, 
they  cannot  be  greater  than  one  nor  less  than  zero.  Moreover,  Q  is 
equal  to  one  when  T'  equals  zero  and  continuously  decreases  as  T' 
increases;  R  is  equal  to  one  when  W  is  zero;  and  for  any  constant 
value  of  T',  R  is  a  continuously  decreasing  function  of  ¥.  Finally, 

S  is  equal  to  one  when  D  is  zero  and,  for  any  particular  pair  of 
values  of  T'  and  ¥,  S  is  a  continuously  decreasing  function  of  D. 

The  point  at  which  the  equiprobability  surface  intersects  the  T'  axis 
Is  determined  by  the  value  of  T'  for  which  Q  =  Pe,  since  R  and  S 
are  both  equal  to  one.  Similarly,  the  surface  intersects  the 
¥  and  D  axes  at  values  of  ¥  and  D  which  make  R  and  S,  respectively, 
equal  to  Pe.  A  constant -temperature  contour  passing  through  the 
point  0X  on  the  equiprobability  surface  would  be  defined  by  the  con¬ 
dition  T’  =  Tx’,  from  which  it  follows  that  Q  =  Qlf  and  hence 
RS  =  Pe/Qi*  These  values  are  depicted  in  figure  2. 

The  preceding  theory  concerning  exceedance  probability  must  be 
modified  slightly  to  make  it  applicable  to  the  particular  problem  under 
consideration,  in  which  the  probabilities  are  to  apply  to  the  most 
severe  icing  conditions  in  an  icing  encounter.  In  most  cases,  the  most 
severe  icing  condition  is  represented  by  the  maximum  value  of  liquid- 
water  content  and  the  corresponding  (simultaneously  observed)  value  of 
drop  diameter.  For  the  portion  of  the  equiprobability  surface  which 
represents  low  values  of  liquid-water  concentration  a.nfl  large  values  of 
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drop  diameter,  on  the  other  hand,  the  critical  condition  is  represented 
by  the  maximum  value  of  drop  size  and  the  corresponding  value  of  liquid- 
water  content.  To  take  account  of  the  first  of  these  possible  combina¬ 
tions,  the  function  R(W,T')  was  determined  from  the  greatest  values  of 
liquid -water  content  for  each  icing  encounter  and  the  function  S(D,W,TT ) 
was  obtained  from  the  corresponding  values  of  drop  diameter;  and  to  take 
account  of  the  second  possible  combinations,  another  function  S'(D,T’) 
was  determined  from  the  greatest  values  of  drop  diameter  and  a  corre¬ 
sponding  function  R'(W,D,T')  was  obtained  from  the  corresponding  values 
Of  liquid-water  content.  The  functions  R(W,T')  and  S(D,Tr)  were  then 
used  to  determine  the  main  portion  of  the  equiprobability  surfaces, 
while  the  functions  R'(W,D,T’)  and  S*(D,T')  were  used  for  the  portion 
representing  large  values  of  D  and  small  values  of  W. 

Construction  of  the  equiprobability  surfaces.-  The  procedure  used 
in  the  construction  of  equiprobability  surfaces  from  a  particular  set  of 
data  was  essentially  that  of  determining  graphically  the  relations  des¬ 
cribing  the  functions  Q(T'),  R(W,T')>  S(D,W,T')>  R'(W,D,T’),  and  S'(D,T*). 
This  procedure  is  illustrated  in  appendix  A. 

The  mathematical  functions  used  to  represent  the  various  distribu¬ 
tions  were  chosen  mainly  on  the  basis  of  the  empirical  criterion  of 
goodness  of  fit,  although  theoretical  considerations  also  were  given 
considerable  weight,  particularly  in  the  selection  of  a  function  to 
represent  the  distribution  of  maximum  values  in  those  cases  in  which 
each  item  in  the  distribution  is  the  maximum  of  a  group  of  observed 
values.  The  function  chosen  for  representing  such  distributions  is  one 
originated  by  Dr.  E.  J.  Gumbel  (reference  11)  to  describe  a  distribution 
of  extreme  values.  This  function  will  be  referred  to  hereinafter  as 
"Gumbel1 s  distribution." 

Because  of  an  almost  complete  lack  of  observational  data  in  flight 
concerning  icing  conditions  at  very  low  temperatures,  reliance  has  been 
placed  on  the  results  of  laboratory  experiments  on  the  formation  of  ice 
crystals  to  aid  further  in  the  construction  of  the  equiprobability  charts. 
The  experiments  relied  upon  (summarized  in  reference  12)  indicate  that  a 
critical  temperature  exists  in  the  neighborhood  of  -^0°  F,  below  which 
ice  crystals  are  formed  in  very  great  numbers  in  atmospheric  air  when¬ 
ever  conditions  of  saturation  with  respect  to  liquid  water  occur.  On 
the  basis  of  this  fact,  it  is  inferred  that  clouds  at  temperatures 
below  -40°  F  are  composed  either  entirely  or  almost  entirely  of  ice 
crystals;  hence  the  liquid-water  content  is  assumed  to  be  zero  at  all 
values  of  T'  greater  than  72°  F. 

The  methods  outlined  in  the  preceding  paragraphs  have  been  used  to 
prepare  constant-temperature  contours  defining  the  equiprobability  sur¬ 
faces  for  values  of  Pe  equal  to  0.1,  0.01,  and  0.001  for  cases  1 
through  5  (table  I).  These  curves  are  presented  in  figures  3  through  7. 
Details  of  the  probability  calculations  for  cases  2  and  5  are  given  as 
examples  in  appendix  A. 
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The  equiprobability  surfaces  defined  by  the  contours  presented  in 
figures  3  through  7  are  constructed  only  for  three  exceedance  probabil¬ 
ities,  but  it  may  be  desired  to  ascertain  a  point  which  lies  between 
the  surfaces  defined  for  the  three  different  exceedance  probabilities. 

If  such  a  point  is  desired,  it  can  be  identified  in  T’-D-W  space  by 
logarithmic  interpolation  between  the  equiprobability  surfaces,  since 
the  probability  functions  decrease  approximately  exponentially  with 
increasing  values  of  T’,  W,  and  D. 

The  inclusion  of  horizontal  extent  as  an  additional  variable. - 
Experience  has  shown  that  the  more  localized  a  cloud  formation,  the 
higher  will  be  its  liquid-water  content.  The  results  of  the  foregoing 
analysis  (fjigs.  3  through  7 ) ,  however,  are  directly  applicable  only  to 
conditions  not  differing  greatly  in  horizontal  extent  from  the  average 
distance  over  which  the  measurements  were  made,  namely,  3  miles  in 
cumulus  clouds  and  10  miles  in  layer  clouds.  The  results  may  be  made 
applicable,  though,  to  other  values  of  horizontal  extent  by  making  use 
of  the  data  contained  in  reference  3  concerning  the  relation  between 
maximum  average  liquid-water  content  and  distance  along  the  flight  path. 

In  order  to  modify  values  of  liquid-water  content  obtained  for  one 
horizontal  extent  to  be  applicable  to  another,  it  has  been  found  neces¬ 
sary  to  make  three  assumptions.  These  assumptions  are  that  (l)  the 
variables,  drop  diameter,  and  temperature  are  assumed  not  to  be  signif¬ 
icantly  related  to  horizontal  extent;  (2)  the  variation  of  average 
liquid-water  content  with  horizontal  extent  is  assumed  to  be  independent 
of  drop  size  and  temperature;  and  (3)  the  results  from  reference  3  are 
assumed  to  be  representative  of  conditions  as  they  would  be  encountered 
in  normal  flight  operations.  This  last  assumption  is  uncertain  because 
of  the  limited  number  of  flights  analyzed  in  reference  3  (11  flights  in 
layer -type  clouds  and  26  in  cumulus  clouds )  and  because  the  flight  paths 
were  frequently  chosen  with  the  aim  of  prolonging  and  maximizing  the 
exposure  to  icing  conditions.  On  the  basis  of  these  assumptions,  equi¬ 
probability  surfaces  applicable  to  other  values  of  horizontal  extent  may 
be  constructed  by  multiplying  all  values  of  liquid-water  content  from 
the  probability  curves  of  figure  3  to  7  Inclusive  by  the  variable  fac¬ 
tor  F  presented  in  figure  8. 

From  figure  8,  it  can  be  seen  that  the  horizontal  extent  selected 
as  being  applicable  to  a  particular  icing  encounter  can  have  a  consid¬ 
erable  effect  on  the  factor  F  used  to  modify  liquid-water-content 
values.  This  effect  is  particularly  pronounced  in  the  case  of  layer 
clouds  where  the  horizontal  extent  can  be  very  large.  For  example,  if 
an  icing  encounter  in  layer  clouds  is  considered  to  extend  over  100  miles 
instead  of  10,  the  factor  to  use  to  modify  liquid-water-content  values 
is  about  0.4  instead  of  1.0.  Consequently,  it  is  necessary  to  consider 
the  appropriateness  for  design  purposes  of  the  standard  values  of  cloud 
horizontal  extent  (3  miles  for  cumulus  and  10  miles  for  layer  clouds). 

For  the  "instantaneous*'  and  "intermittent"  classifications  of  reference  5 
(classes  I  and  II  as  defined  therein),  the  horizontal  extents 
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of  l/2  and  3  miles,  respectively,  as  assigned  in  reference  5,  are  con¬ 
sidered  to  be  reasonable  values  for  design  purposes.  For  the  "contin¬ 
uous"  classification  (class  III),  a  horizontal  extent  of  20  miles 
appears  appropriate  for  design.  Although  it  is  realized  that  contin¬ 
uous  icing  might  well  extend  beyond  20  miles,  reduction  of  the  liquid- 
water  content  by  a  factor  less  than  0.8  (an  approximate  value  of  F 
corresponding  to  20  miles  horizontal  extent )  does  not  appear  justified 
for  design  purposes.  It  should  be  noted  that  these  three  design  values 
of  horizontal  cloud  extent,  shown  in  figure  8  by  broken  lines,  are  sug¬ 
gested  values  only  and  do  not  preclude  the  use  of  the  curves  for  the 
estimation  of  the  factor  F  by  selecting  other  values  of  horizontal 
extent . 


The  Determination  of  the  Probability  of  Exceeding  Any  Specified 
Group  of  Values  of  Liquid-Water  Content  Associated  Simul¬ 
taneously  With  Temperature  and  Drop-Diameter  Values 
Lying  Within  Specified  Ranges 


The  equiprobability  surfaces  discussed  in  the  preceding  section 
provide  information  concerning  the  simultaneous  combinations  of  liquid- 
water  content,  drop  diameter,  and  temperature  having  the  same  probability 
of  occurrence.  Another  useful  expedient  for  the  designer,  however,  would 
be  a  means  for  determining  the  probability  of  encountering  any  icing 
condition  not  included  by  the  values  of  liquid-water  content,  drop  diam¬ 
eter,  and  temperature  specified  as  design  criteria  for  a  particular  com¬ 
ponent  of  a  thermal  system.  The  basis  for  such  a  technique  lies  in 
determining  the  probability  of  equaling  or  exceeding  certain  values  of 
liquid-water  content  specified  for  a  number  of  icing  conditions. 

Probability  of  encountering  icing  conditions  which  exceed  certain 
specified  values.-  If,  for  a  particular  element  of  a  thermal  system  and 
a  chosen  value  of  horizontal  cloud  extent,  combination  values  of  liquid- 
water  content,  drop  diameter,  and  temperature  are  specified  which  the 
element  can  just  tolerate,  the  values  would  define  a  surface  in 
T'-D-W  space.  Such  a  surface  can  be  viewed  as  representing  the  locus 
of  liquid-water-content  values  to  which  the  element  is  critical  over 
specified  ranges  of  temperatures  and  drop  diameters.  Individual  values 
on  such  a  surface  would  be  obtained  through  knowledge  of  the  area,  rate, 
and  distribution  of  water-drop  impingement  for  the  particular  element. 

For  purposes  of  explanation,  figure  9  is  presented  to  show  the 
locus  of  critical  values  of  liquid-water  content  for  an  assumed  element 
of  a  thermal  system.  The  locus  of  liquid-water-content  values  will  be 
termed  a  "marginal"  surface  to  indicate  that  values  of  W  lying  above 
the  surface  cannot  be  tolerated  by  the  element  of  the  thermal  system 
under  consideration.  In  general,  the  shape  of  the  marginal  surface  in 
the  W-D  plane  for  the  case  of  a  wing  element  having  a  given  heat  out¬ 
put  may  be  assumed  to  have  a  shape  resembling  a  hyperbola.  A  hyperbolic 
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shape  is  appropriate  because  for  drop  diameters  less  than  the  critical 
diameter4  for  the  element  under  consideration,  an  infinitely  large 
amount  of  liquid  water  is  tolerable;  and,  as  drop  diameter  progressively 
increases,  decreasing  amounts  of  liquid  water  are  tolerable.  In  any 
W-T’  plane,  however,  W  will  be  assumed  to  vary  from  the  values 
of  W  at  T'  =  0  to  ¥  =  0  at  some  value  of  Tx',  say  Tcr',  where 
the  heating  capacity  of  the  element  can  bring  the  surface  temperature 
of  the  element  just  up  to  freezing  temperature  (a  clear-air  condition). 

For  convenience,  the  base  of  the  marginal  surface  (the  T'-D  plane 
at  W  =  0)  is  divided  into  n  rectangles  of  dimension  AD  and  AT'  such 
that  variations  are  small  in  values  of  liquid-water  content  included  by 
the  projection  of  each  rectangle  on  the  marginal  surface.  The  value  of 
liquid-water  content  considered  to  be  representative  for  each  rectangle 
is  designated  as  Wp  in  figure  9»  The  values  of  T',  T,  and  D  corre¬ 
sponding  to  Wp  shall  be  called  Tp',  Tp,  and  Dp,  respectively.  For 
every  set  of  Wp-Dp-Tp'  values,  there  exists  a  probability  function  APp 
which  expresses  the  probability  that,  for  a  random  icing  encounter,  the 
liquid-water  content  will  exceed  Wp  and  the  drop  diameter  and  tempera¬ 
ture  depression  will  be  within  the  ranges  AD  and  AT',  respectively. 

The  n  values  of  APp  corresponding  to  the  n  rectangles  in  the 
T'-D  plane  represent  the  probabilities  of  the  occurrence  of  an  event 
(the  overloading  of  the  anti-icing  element)  in  n  different,  mutually 
exclusive  ways.  The  probability  that  the  event  will  occur  is  the  sum 
of  the  probabilities  that  it  will  occur  in  any  number  of  different  ways 
provided  the  different  possibilities  are  mutually  exclusive.  (See 
reference  10.)  Hence,  the  probability  of  overloading  the  given  anti¬ 
icing  element  is  approximately  equal  to  the  sum  of  the  partial  probabil¬ 
ities  for  the  n  rectangles  according  to  the  equation 

P  =  |  APp  (3) 

In  order  to  evaluate  APp  for  each  of  the  n  rectangles,  it 
would  appear  that  a  separate  determination  of  the  relationship 
between  W  and  APi  for  each  rectangle  would  be  necessary.  This  pro¬ 
cedure  is  not  practicable,  however,  because  if  the  number  of  rectangles 
is  chosen  large  enough  to  make  variations  in  W  small  for  each  rec¬ 
tangle,  then  the  amount  of  data  available  to  define  APp  as  a  function 
of  W  is  too  small  to  yield  statistically  reliable  results  for  the 
individual  rectangles.  Moreover,  the  presentation  of  the  results  would 
require  separate  curves  or  tables  for  each  rectangle.  It  is  possible 
to  overcome  these  difficulties  by  utilizing  the  fact  that  the 

4The  critical  diameter  may  be  defined  as  the  smallest  diameter  for  which 
drops  will  impinge  on  the  element.  For  some  airfoils,  the  critical 
size  can  be  estimated  by  a  technique  outlined  in  NACA  TU  2kj6  (refer- 
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distribution  of  drop  diameter  is  approximately  independent  of  tempera¬ 
ture  and  liquid-water  content.  Although  there  is  a  slight  tendency 
toward  higher  frequencies  of  large  values  of  drop  diameter  when  the  w 

liquid-water  content  is  low,  this  effect  is  too  small  to  have  an  impor¬ 
tant  influence  in  the  determination  of  the  total  marginal  probability. 

By  maki ng  the  assumption  that  the  distribution  of  drop  diameter  is 
independent  of  temperature  and  liquid-water  content,  it  is  possible  to 
determine  values  of  APp  for  different  combinations  of  temperature 
depression  below  freezing,  liquid-water  content,  drop  diameter,  and 
horizontal  extent. 

When,  for  any  given  cloud  horizontal  extent,  drop-diameter  distri¬ 
bution  is  assumed  to  be  independent  of  temperature  and  liquid  water,  the 
probability  equation 

APi  =  P^  PAT  PWi(AT)  (4) 

can  be  written  from  analogy  with  equation  (l).  The  interpretation  of 
equation  (4)  is  that  the  partial  probability  APp  applicable  to  a  - 

chosen  point  on  the  marginal  surface  is  the  product  of  three  probabil¬ 
ities:  the  probability,,  P^p,  that  the  drop  diameter  will  be  within  the 
interval  AD;  the  probability,  P^rp,  that  the  temperature  will  lie  within  w 

the  interval  AT;  and  the  conditioned  probability,  Pjtf.(AT),  that  the 
liquid-water  content  will  be  equal  to  or  greater  than  Wj_  under  the 
condition  that  the  temperature  lies  within  the  interval  AT. 

To  evaluate  equation  (4),  values  of  P^j)  were  obtained  from  the 
frequency  distribution  of  all  observations  of  drop  diameter;  values 
of  Pat  were  obtained  from  the  frequency  distribution  of  the  tempera¬ 
ture  of  the  most  severe5  icing  observations  per  icing  encounter;  and 
values  of  P^(AT)  were  obtained  from  the  frequency  distributions  of 

maximum  liquid-water  content  per  encounter  for  the  various  temperature 
intervals.  Because  of  the  assumption  that  drop  size  is  independent  of 
temperature  and  liquid -water  content,  it  was  not  possible  to  provide 
exactly  for  the  case  in  which  the  maximum  observed  drop  diameter  deter¬ 
mines  the  most  severe  icing  condition  in  an  icing  encounter.  An  approx¬ 
imate  allowance  was  made  for  this  factor,  however,  by  substituting  the 
frequency  distribution  of  the  greatest  value  of  drop  diameter  per 
encounter  over  a  limited  range  of  large  drop  diameters. 

Equation  (4)  was  evaluated  for  different  horizontal  extents,  and  to 
do  this  the  frequency  distribution  for  maximum  liquid-water  content  per 
encounter  was  altered  according  to  a  relation  from  reference  3  between 
maximum  average  liquid-water  content  per  encounter  and  the  distance 
along  the  flight  path. 

* 


5For  the  purposes  of  this  report,  the  most  severe  icing  condition  is  con- 
sidered  as  the  one  having  the  largest  liquid-water  content. _ 
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The  partial  probability  charts.-  In  evaluating  equation  (4)  for 
any  one  value  of  cloud  horizontal  extent,  numerous  values  of  APj_  are 
possible  from  the  meteorological  data  for  different  combinations 
of  P/£),  Pat j  and  P^(AT).  The  values  of  APp,  however,  can  be  conven¬ 
iently  summarized  in  the  form  of  partial  probability  charts.  Such 
charts  are  presented  in  figure  10  for  cases  1  through  5  of  table  I. 

As  is  evident  from  figure  10,  seven  intervals  of  drop  diameter  and 
six  intervals  of  temperature  were  used  in  the  construction  of  the 
charts.  These  intervals  were  chosen  so  as  to  avoid  large  variations 
of  W  encompassed  by  any  particular  AD-AT’  rectangle  projected  on 
the  marginal  surface.  In  all  conceivable  combinations,  the  temperature 
and  drop-diameter  intervals  provide  42  different  values  of  APp  which 
must  be  summed  (in  accordance  with  equation  (3))  to  obtain  the  probabil¬ 
ity,  P,  of  overloading  the  particular  element. 

To  obtain  a  specific  value  of  APp  from  a  chart  of  figure  10,  the 
following  procedure  is  used:  Choose  a  value  of  liquid -water  content  Wp 
and  follow  a  vertical  line  until  it  intersects  with  the  desired 
horizontal-extent  line;  from  this  intersection,  follow  horizontally  until 
an  intersection  with  the  appropriate  temperature -interval  curve  is 
obtained;  from  the  point  of  this  intersection,  follow  a  vertical  line 
downward  to  obtain  an  intersection  with  a  chosen  drop-diameter  curve; 
a  horizontal  line  through  the  latter  intersection  yields  the  desired 
value  of  APp  on  the  ordinate  scale.  For  ease  in  observing  this  pro¬ 
cedure,  each  chart  in  figure  10  has  a  system  of  arrows  leading  from 
some  arbitrary  value  of  Wp  through  the  horizontal-extent,  temperature, 
and  drop-diameter  curves  to  a  corresponding  value  Of  APp. 

It  should  be  noted  that  partial  probabilities  presented  in  figure  10 
can  be  determined,  for  values  of  horizontal  extent  not  specifically 
represented  in  the  charts,  by  interpolation  between  the  horizontal-extent 
curves.  It  is  not  permissible,  however,  to  interpolate  between  the 
curves  representing  temperature  and  drop  diameter  since  the  partial  prob¬ 
abilities  are  dependent  upon  the  size  of  the  intervals  for  which  the 
temperature  and  diameter  curves  are  constructed.  The  approximate  median 
values  of  temperature  and  drop  diameter,  which  correspond  to  the  intervals 
of  Tp  and  Dp  used  In  the  charts  of  figure  10,  are  listed  in  table  II 
for  reference.  Median  values  are  used  because  the  median  has  the  advan¬ 
tage  that  its  value  is  not  greatly  influenced  by  the  chance  variations 
of  extreme  items.  A  detailed  description  of  the  techniques  used  in 
constructing  the  partial  probability  charts,  with  application  to  selected 
cases,  Is  presented  in  appendix  B. 
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DISCUSSION 

Reliability  of  the  Equiprobability  and  the 
Partial  Probability  Charts 

The  reliability  of  the  probability  analysis  reported  herein  is 
affected  chiefly  by  four  factors;  namely,  (l)  errors  of  measurement, 
(2)  amount  of  data  available,  (3)  representativeness  of  the  data,  and 
(4)  limitations  with  respect  to  climate  and  altitude.  The  possible 
effect  of  these  factors  is  discussed  in  the  following  paragraphs. 


Errors  of  measurement.-  As  noted  in  reference  l4,  the  measurement 
of  liquid-water  content  and  drop  diameter  by  the  multicylinder  method 
is  subject  to  error.  The  following  factors  may  be  listed  as  possible 
sources  of  error:  (l)  runoff,  (2)  bounce-off,  (3)  evaporation, 

(4)  variations  in  drop  concentration  in  the  neighborhood  of  the  cylinders 

due  to  disturbing  influences  of  the  airplane,  (5)  collection  of  snow 

flakes  as  well  as  liquid  drops  in  mixed  clouds,  (6)  inaccuracies  due  to 

the  theoretical  drop-size  distributions  differing  from  those  assumed  in 

the  calculations,  (7)  failure  to  fulfill,  with  the  apparatus  used,  the 

theoretical  assumption  of  two-dimensional  flow,  (8)  irregularities  in  ^ 

the  surfaces  of  the  ice-covered  cylinders,  (9)  errors  in  weighing, 

(10)  irregular  breaking  of  the  ice  when  separating  cylinders,  (ll)  errors 
in  measuring  the  final  diameter  of  the  ice-covered  cylinders  or  in  cal¬ 
culating  the  average  diameter  from  an  assumed  value  of  ice  density, 

(12)  errors  in  the  measurement  of  airspeed  and  exposure  time,  and 

(13)  errors  in  matching  the  experimental  data  with  the  theoretical 
collection-efficiency  curves . 

A  detailed  discussion  of  the  possible  effects  of  each  of  these 
factors  on  the  accuracy  of  the  probability  analysis  presented  herein  is 
beyond  the  scope  of  this  report.  The  problem  of  runoff,  though,  is 
worthy  of  discussion  since  this  factor  imposes  an  upper  limit  to  the 
range  of  liquid-water-content  values  that  can  be  measured.  The  results 
of  heat -transfer  calculations  reported  in  reference  15  indicate  that  the 
maximum  value  of  liquid-water  content  that  can  be  measured  with  rotating 
cylinders  as  they  are  normally  used  varies  approximately  linearly  with 
temperature  depression  below  freezing,  reaching  about  2  grams  per  cubic 
meter  at  5°  F  for  an  airspeed  of  about  200  miles  per  hour.  An  examina¬ 
tion  of  the  data  used  in  the  present  analysis  indicates  that  not  more  than 
about  5  percent  of  the  icing  encounters  in  any  of  the  cases  analyzed  con¬ 
tained  observed  values  of  liquid-water  content  which  were  close  enough  to 
the  theoretical  maximums  to  be  regarded  as  likely  to  be  influenced  by 
runoff.  In  order  to  estimate  the  effect  of  a  small  percentage  of  errors 
of  this  type  on  the  results  of  the  probability  analysis,  it  is  necessary 
to  consider  the  manner  in  which  the  liquid-water-content  distribution  a 

functions  were  determined  from  the  experimental  data.  In  this  process, 
the  cumulative  frequency  distribution  of  maximum  liquid-water  content 
per  icing  encounter  was  plotted  on  Gumbel's  distribution  paper,  and  the 
straight  line  of  best  fit  was  drawn  to  represent  the  entire  distribution; 
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thus  errors  in  the  upper  5  percent  of  the  distribution  had  only  a  minor 
effect  in  the  final  selection  of  the  line  used  to  represent  the  distri¬ 
bution  function.  Moreover,  it  was  noted  that  in  nearly  al 1  cases  the 
entire  distribution  followed  the  straight  line  with  minor  and  apparently 
random  variations  and  without  any  consistent  or  marked  tendency  for  large 
negative  departures  at  the  upper  end  of  the  distribution.  It  may  be 
inferred,  therefore,  that  the  limitation  on  the  measurement  of  liquid- 
water  content  discussed  in  reference  15  does  not  have  an  important  influ¬ 
ence  on  the  results  of  this  analysis. 

In  view  of  the  number  and  complexity  of  the  possible  sources  of 
error,  it  is  not  possible  at  this  time  to  make  a  reliable  estimate  of  the 
total  accuracy  of  the  measurement  of  liquid-water  content.  In  a  statis¬ 
tical.  analysis,  however,  reasonably  reliable  results  may  be  obtained  as 
long  as  the  errors  are  small  compared  to  the  real  variations  of  the 
quantity  measured  and  are  distributed  approximately  at  random.  It  is 
believed  that  these  conditions  are  approximately  fulfilled  for  the  values 
of  liquid-water  content  measured  by  the  multicylinder  technique. 

Multicylinder  measurements  of  drop  diameter,  on  the  other  hand,  while 
quite  accurate  and  reliable  for  small  drops,  become  increasingly  inaccurate 
as  the  drop  size  increases.  Moreover,  the  errors  are  not  normally  dis¬ 
tributed,  since  large  positive  errors  are  more  probable  than  large  nega¬ 
tive  errors,  especially  at  large  values  of  drop  diameter  (see  reference  l4) 
The  result  of  these  errors  is  an  increase  in  the  dispersion  of  the  fre¬ 
quency  distribution  and  an  exaggeration  of  the  probability  of  occurrence  of 
large  values  of  drop  diameter.  This  effect  is  probably  not  important  for 
observations  in  the  eastern  United  States  because  in  this  case  the  drop 
size  rarely  exceeds  the  range  of  reliable  measurement.  In  the  Pacific 
coast  area,  on  the  other  hand,  where  large  drops  occur  with  much  greater 
frequency,  it  is  quite  likely  that  the  extreme  values  of  drop  diameter  are 
influenced  considerably  by  this  factor. 

Amount  of  available  data.-  Some  concept  of  the  accuracy  of  the  proba¬ 
bility  analyses  for  the  various  cases  may  be  obtained  by  reference  to  the 
number  of  icing  encounters  for  each  geographical  area  and  each  cloud-type 
classification.  These  data  are  presented  in  table  I.  It  Is  noted  that 
the  greatest  amount  of  data  for  any  one  case  is  110  encounters  in  layer 
clouds  in  the  eastern  United  States.  A  sample  of  this  size  is  sufficient 
to  establish  the  probabilities  with  considerable  confidence  at  the  0.01 
probability  level,  and  the  extrapolation  to  a  probability  of  0.001  is  not 
likely  to  introduce  serious  errors.  In  the  remaining  cases,  the  sample 
size  ranges  from  25  to  Ml-  encounters.  In  these  cases,  the  extrapolation 
to  a  probability  of  0.001  may  introduce  considerable  uncertainty. 

Representativeness  of  the  data.-  Another  factor,  which  probably  has 
as  great  an  influence  on  the  reliability  and  applicability  of  the  analy¬ 
sis  as  either  errors  of  measurement  or  sample  size,  is  the  systematic  bias 
Introduced  by  the  fact  that  all  the  test  data  were  obtained  during  flights 
in  which  icing  conditions  were  deliberately  sought  instead  of  being  acci¬ 
dentally  encountered.  It  is  difficult  to  estimate  the  effect  of  this 
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factor  since  measurements  of  liquid-water  content  and  cloud-drop  size 
encountered  during  ordinary  transport  operations  are  not  yet  available. 

It  would  appear  that  the  principal  difference  between  transport  experi¬ 
ence  and  icing  research  flights  would  be  the  greater  frequency  of  encoun¬ 
ters  with  icing  conditions  in  the  latter.  This  difference  alone,  however, 
would  not  affect  the  validity  of  the  analysis  since  the  icing  encounter  is 
taken  as  the  unit  of  experience  upon  which  the  analysis  is  based.  The 
question'  to  be  considered  is  whether  the  icing  encounters  which  occurred 
during  the  research  flights  represent  an  unbiased  sample  of  the  encoun¬ 
ters  which  would  occur  during  normal  flight  operations. 

The  differences  in  severity  and  extent  between  icing  conditions 
experienced  during  research  flights  and  the  conditions  which  would  have 
been  encountered  during  an  ordinary  flight  in  the  same  general  area  at 
the  same  time  are  of  two  types.  The  first  type  occurred  in  cases  in 
which  no  icing  would  have  been  encountered  on  a  normal  flight  through 
the  area  but,  due  to  the  existence  of  clouds  in  thin  layers  or  localized 
areas it  was  possible  to  experience  icing  with  the  test  airplane. 
Encounters  of  this  type  would  be  expected  to  be  less  severe  and  exten¬ 
sive,  on  the  average,  than  typical  encounters  during  normal  operations. 
Differences  of  the  second  type  occurred  in  cases  in  which  icing  condi¬ 
tions  would  have  been  encountered  on  a  normal  flight  but,  due  to  a  con¬ 
scious  effort  on  the  part  of  the  pilot,  the  conditions  experienced  by 
the  research  airplane  were  more  severe  and  extensive.  It  is  evident, 
therefore,  that  the  special  procedures  followed  during  research  flights 
have  resulted  in  the  inclusion  in  the  test  results  of  data  from  icing 
encounters  which  were  both  more  and  less  severe  than  those  which  would 
have  been  experienced  under  normal  operating  conditions.  It  is  not 
possible  to  state  with  certainty  what  is  the  net  effect  of  these  two 
types  of  differences,  but  it  is  believed  that  differences  of  the  second 
type  predominate,  with  the  result  that  the  data  are  somewhat,  though  not 
greatly,  biased  by  the  inclusion  of  an  abnormally  high  frequency  of 
severe  icing  conditions. 

Effect  of  climate  and  altitude. -  Another  factor  having  an  important 
bearing  on  the  applicability  of  the  results  is  the  effect  of  climate  and 
flight  altitude  on  the  temperature  of  icing  conditions.  It  is  likely 
that  the  distribution  of  temperatures  of  icing  conditions  encountered  by 
the  test  airplanes  closely  resembles  that  which  would  be  found  in  normal 
operations  of  unpressurized  airplanes  during  winter  and  spring  in  the 
areas  studied.  Different  temperature  distributions  would  be  expected, 
however,  with  high-altitude  airplanes  or  with  conventional  airplanes  in 
other  climateB  and  seasons.  For  example,  if  an  airplane  is  expected  to 
cruise  at  altitudes  of  from  18,000  to  28,000  feet,  where  the  temperature 
Is  normally  In  the  range  between  0  and  -40°  F,  a  much  higher  frequency 
of  icing  at  low  temperatures  is  to  be  expected  and  the  results  of  an 
analysis  of  icing  conditions  occurring  mostly  at  higher  temperatures 
would  be  applicable  only  to  a  limited  extent. 
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Comparison  of  Results  of  the  Equiprobab ill ty 
Analysis  With  Results  of  MCA  TN  1855 


Certain  values  of  liquid-water  content,  drop  diameter,  and  temper¬ 
ature  were  presented  in  MCA  TN  1855  (reference  5)  for  consideration  in 
the  design  of  anti-icing  equipment.  The  combinations  of  meteorological 
variables  recommended,  however,  were  lacking  in  one  respect,  namely,  the 
probability  of  occurrence  of  these  particular  combinations.  Fortunately, 
with  the  aid  of  the  equiprobability  curves,  some  indication  of  these 
probabilities  can  be  obtained  when  corresponding  meteorological  condi¬ 
tions  are  compared. 

Four  classes  of  icing  conditions  are  presented  in  MCA  TN  1855* 

I,  instantaneous ;  II,  intermittent;  III,  continuous;  and  IV,  freezing 
rain.  All  classes,  except  IV,  are  subdivided  into  two  types  of  icing 
conditions,  maximum  and  normal,  and  are  confined  to  icing  conditions 
associated  with  definite  cloud  formations  of  characteristic  horizontal 
extents.  Hence,  classes  I,  II,  and  III  are  ideal  for  comparison  with 
the  results  which  can  be  obtained  from  the  equiprobability  charts 
(figs.  3  to  7)  and  the  relations  for  horizontal  extent  (fig.  8). 

Two  groups  of  icing  conditions  are  chosen  as  a  basis  for  comparison 
between  the  analysis  of  this  report  and  NACA  TN  1855.  The  first  group 
consists  of  icing  conditions  selected  from  the  probability  analysis  for 
the  case  of  layer-type  clouds  with  an  exceedance  probability 
of  Pe  =  0.001.  The  second  group  consists  of  icing  conditions  for  the 
case  of  Pacific  coast  cumulus,  also  with  an  exceedance  probability 
of  0.001.  These  two  groups  are  compared,  respectively,  with  classes  III-M 
and  II-M  from  NACA  TN  I855.  The  comparisons  are  made  by  determining 
values  of  liquid-water  content  from  figure  4  and  figures  3,  5,  and  7  for 
the  various  combinations  of  temperature  and  drop  diameter  given  in 
table  I  of  MCA  TN  1855  for  class  II-M  and  III-M  conditions,  respectively. 
The  three  values  of  liquid-water  content  obtained  from  figures  3,  5, 
and  7  corresponding  to  class  III-M  were  averaged  by  weighting  the  values 
according  to  the  horizontal  extent  of  the  region  to  which  they  apply. 

The  three  regions  (fig.  1)  have  area  ratios  of  about  1/8,  2/8,  and  5/8 
for  the  Pacific  coast,  plateau,  and  eastern  regions,  respectively;  and 
these  area  ratios  were  used  to  weight  the  values  of  liquid-water  content 
from  the  equiprobability  charts  so  that  the  resulting  average  values 
would  be  on  a  comparable  basis  with  the  values  expressed  in  MCA  TN  1855. 
The  weighted  average  values  of  liquid-water  content  were  corrected  by 
figure  8  to  make  them  applicable  to  a  horizontal  extent  of  20  miles,  a 
figure  reasonably  applicable  to  the  layer-cloud  data  of  NACA  TN  1855. 

For  class  II-M  conditions,  no  correction  for  horizontal  extent  was 
required  because  the  values  in  MCA  TN  1855  and  figure  k  both  apply  to 
3  miles  horizontal-cloud  extent. 

The  results  of  the  comparisons  between  NACA  TN  1855  and  the  proba¬ 
bility  analysis  are  shown  in  tables  III  and  IV  for  classes  II-M  and  III-M, 
respectively.  In  general,  the  liquid-water-content  values  agree  quite 
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closely,  except  at  temperatures  of  -4°  F  and  below.  In  this  low  tempera¬ 
ture  range,  the  values  of  liquid-water  content  derived  from  the  equiprob- 
ability  surfaces  are  considerably  lower  than  values  given  in  NACA  TN  1855* 
The  reason  is  that  the  probability  analysis  entails  some  extrapolation  in 
this  temperature  range,  and  so  may  not  provide  exactly  the  correct  value 
of  liquid-water  content.  On  the  other  hand,  the  values  listed  in 
NACA  TN  1855  are  not  restricted  to  actual  measurements,  and  therefore  may 
not  be  of  precise  magnitude.  For  design  purposes,  however,  the  values 
given  by  the  probability  analysis  should  be  of  proper  magnitude.  In  this 
regard,  it  should  be  noted  that  the  majority  of  the  data  utilized  in  the 
probability  analysis  was  taken  at  comparatively  low  altitude  (13,000  feet), 
whereas  the  temperature  range  between  -4°  F  and  -40°  F  represents  consider¬ 
ably  higher  altitudes  (l8,000  to  28,000  feet).  This  difference  in  altitude 
could  have  some  bearing  on  the  accuracy  of  the  probability  data,  particu¬ 
larly  at  very  high  altitudes. 

By  a  procedure  similar  to  that  used  in  the  comparisons  for  the  two 
example  cases,  the  correspondence  between  classes  of  icing  conditions 
represented  by  the  equiprobability  charts  and  other  classes  of  conditions 
presented  in  NACA  TN  1855  could  be  ascertained  approximately.  The  corre¬ 
spondence  for  these  classes  of  conditions,  and  also  for  the  two  example 
cases,  are  shown  in  table  V. 

An  inspection  of  table  V  reveals  a  consistency  between  the. classes 
of  design  values  recommended  in  NACA  TN  1855  and  specific  values  of 
exceedance  probabilities.  Only  the  instantaneous  maximum  icing  condi¬ 
tion  (Pe  =  0.0001)  appears  to  be  incongruous  with  the  maximum  conditions 
in  classes  II-M  and  III-M  (Pe  =  0.001).  This  apparent  Incongruity  is 
directly  attributable  to  the  fact  that  the  instantaneous  maximum  condi¬ 
tion  presented  in  NACA  TN  1855  was  calculated  for  tall  tropical  cumulus 
clouds.  Such  extremely  severe  icing  conditions  would,  of  course,  be 
exceeded  infrequently  as  is  borne  out  by  the  probability  analysis. 

Another  point  of  interest  in  regard  to  table  V  is  the  general  order  of 
magnitude  of  the  exceedance  probabilities  which  are  found  to  apply  to 
the  various  classes  of  icing  conditions  specified  in  NACA  TN  1855*  For 
example,  any  one  of  the  icing  conditions  listed  under  class  II-M  (inter¬ 
mittent  maximum)  or  class  III-M  (continuous  maximum)  will  be  exceeded  once 
in  only  about  a  thousand,  encounters  with  icing  conditions.  It  is  believed 
that  the  icing  conditions  specified  by  classes  II-M  and  III-M  do  not 
impose  too  severe  design  requirements  and,  therefore,  if  the  ice- 
prevention  system  will  cope  with  these  conditions,  the  results  of  the 
probability  analysis  indicate  that  satisfactory  ice  protection  would  be 
provided  for  the  vast  majority  of  icing  encounters. 

USE  OF  THE  EQUIPROBABILITY  AND  PARTIAL  PROBABILITY 
CHARTS  AS  APPLIED  TO  THE  DESIGN  OF  THERMAL  ICE- 
PREVENTI0N  EQUIPMENT 

Heretofore,  no  Information  has  been  available  regarding  the  prob¬ 
ability  of  encountering  icing  conditions  which  exceed  certain  specified 
values,  so  that  some  discussion  of  probability  as  applied  to  the  design 
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of  thermal  anti-icing  equipment  seems  worthwhile.  In  general,  two  design 
problems  exist:  one  problem  is  to  establish  consistent  simultaneous  com¬ 
binations  of  liquid-water  content,  drop  diameter,  and  temperature  which 
can  occur  in  nature,  and  the  other  problem  is  to  determine  how  effectively 
a  given  thermal  element  of  an  anti -icing  system  will  cope  with  the  entire 
array  of  icing  conditions  existing  in  nature.  It  is  the  purpose  of  the 
equiprobability  charts  to  aid  in  the  solution  of  the  first  problem,  while 
the  partial  probability  charts  are  intended  to  aid  in  the  solution  of  the 
second  problem. 


The  Equiprobability  Charts 

The  equiprobability  charts  essentially  are  a  means  to  determine 
quickly  consistent  sets  of  liquid-water  content,  temperature,  and  drop- 
diameter  values  which  have  the  same  probability  of  being  exceeded  during 
icing  encounters.  Sets  of  such  values  obtained  from  the  equiprobability 
charts  can  be  used  as  criteria  for  the  design  of  thermal  anti-icing 
equipment.  No  need  to  use  these  charts  exists  if  it  is  found  that 
table  V  will  define  a  class  of  conditions  from  NACA  TN  1855  suitable  for 
the  needs  of  the  particular  design  problem.  However,  should  it  be 
desired  to  design,  say,  a  jet-engine  inlet  duct  for  an  exceedance  prob¬ 
ability  of  0.001  (only  one  chance  in  a  thousand  of  exceeding  a  specified 
amount  of  liquid  water),  table  V  would  not  be  of  help  because  correspond¬ 
ing  icing  conditions  from  NACA  TN  1855  are  not  defined  for  this  partic¬ 
ular  case.  Hence,  reliance  would  have  to  be  placed  on  the  equiprobability 
charts  in  order  to  establish  consistent  sets  of  values  for  a  chosen 
value  of  Pe . 

General  procedure  for  using  the  equiprobability  charts.-  The  follow¬ 
ing  procedure  for  the  use  of  the  equiprobability  charts  is suggested : 

1.  Determine  for  what  class  of  icing  condition  (I,  II,  or  III  of 
NACA  TN  1855)  a  design  is  to  be  made  by  considering  the  nature  of  the 
component  to  be  protected.  (See  table  Y. ) 

2.  Establish  the  Severity  of  icing  condition  within  the  class  by 
selecting  the  exceedance  probability  for  which  the  component  is  to  be 
designed.  (See  table  V.) 

3*  With  the  aid  of  table  V  and  steps  1  and  2,  determine  what  equi- 
probability  chart  to  use  in  defining  values  of  W  consistent  with  dif¬ 
ferent  sets  of  T  and  D.  The  selection  of  charts  can  be  made  by  noting 
the  type  of  cloud  formation  applicable  to  the  class  (i.e.,  I,  II,  or  III) 
and  severity  of  icing  condition  established  in  steps  1  and  2.  For  cumulus 
clouds,  it  is  suggested  that  the  charts  for  the  Pacific  coast  area  be 
used  for  design  purposes  rather  than  the  charts  for  the  plateau  area. 

The  reason  is  that,  until  more  data  are  available  for  the  eastern 
United  States,  Pacific  coast  cumulus -cloud  data  should  provide  liquid- 
water-content  values  that  are  conservative  for  the  remainder  of  the 
United  States .  For  layer  clouds ,  however,  it  is  suggested  that  the 
charts  for  all  three  geographical  areas  be  used  for  design  purposes  by 
using  a  weighted  average  of  liquid-water-content  values  from  the  three 
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charts.  The  justification  for  using  a  weighted  average  is  that  layer 
clouds  can  he  regarded  as  being  structurally  continuous  over  exceedingly 
large  geographical  areas;  consequently,  only  average  values  of  the 
parameters  defining  sin  icing  condition  would  seem  to  have  much  signif¬ 
icance  in  the  design  of  an  anti-icing  component  for  layer  clouds  exist¬ 
ing  throughout  the  United  States. 

4.  From  the  appropriate  equiprobability  chart,  choose  values 
of  W  corresponding  to  a  variety  of  sets  of  T  and  D.  In  the  case  of 
layer  clouds,  where  more  than  one  chart  may  be  employed  to  determine 
values  of  W,  a  suitable  weighted  average  for  the  entire  United  States 
can  be  obtained  by  weighting  the  liquid-water-content  values  for  the 
Pacific  coast,  plateau,  and  east  coast  regions,  in  the  ratios  of  1/8,  2/8 
and  5/8,  respectively.  These  weight  ratios,  l/8,  2/8,  and  5/8,  are  pro¬ 
portional  to  the  three  geographical  areas  and  may  be  considered  as 
approximately  accounting  for  the  length  of  time  that  an  airplane  flying 
transcontinentally  would  be  in  each  area. 

5-  Multiply  the  values  of  W  determined  in  step  4  by  a  factor, 
chosen  from  figure  8,  for  including  the  effect  of  horizontal  extent  of 
cloud  formation.  The  factor  chosen  should  lie  within  the  range  of 
horizontal  extent  indicated  in  figure  8  for  the  particular  class  icing 
condition  established  in  step  1;  it  should  correspond  also  to  the  curve 
for  the  exceedance  probability  decided  upon  in  step  2.  Recommended 
values  of  the  factor  correspond  to  values  of  horizontal  extent  identi¬ 
fied  in  figure  8  as  design  values.  These  design  values  are  for  cloud 
horizontal  extents  of  1/2,  3;  and  20  miles  and  correspond  to  icing  con¬ 
dition  classes  I,  II,  and  III,  respectively. 

Example  showing  the  use  of  the  equiprobability  charts.-  For  this 
example,  a  design  of  a  jet-engine  inlet-guide  vane  will  be  assumed  in 
which  it  is  desirable  to  protect  against  those  icing  conditions  having 
a  probability  of  being  exceeded  once  in  a  thousand  icing  encounters 
(Pe  =  0.001).  The  steps  in  determining  sets  of  values  of  W,  T,  and  D 
are: 


1.  From  table  V,  jet-engine  intakes  are  recommended  to  be  designed 
for  class  I  icing  conditions. 

2.  The  exceedance  probability  for  which  the  design  is  to  be  made 
is  Pe  =  0.001. 


3.  In  table  V,  Pacific  coast  cumulus  clouds  are  considered  typical 
for  class  I  icing  conditions.  Hence,  equiprobability  charts  for 
Pacific  coast  cumulus  (fig.  4(c))  are  chosen  for  obtaining  values  of 
liquid-water  content  corresponding  to  various  combinations  of  temperature 
and  drop  diameter. 

4.  Values  of  W  selected  from  figure  4(c)  for  different  arbitrar¬ 
ily  chosen  sets  of  T  and  D  are  presented  in  the  following  table: 
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Instantaneous  Conditions  Corresponding  to 
Pe  =0.001 

[Horizontal  extent  not  considered] 


2.65 

32 

2.39 

1.91 

1.20 

•55 

.  . . 

-20 

1.50 

32 

1.02 

10 

.90 

0 

.49 

-10 

.18 

-20 

•  57 

32 

.49 

10 

•  32 

•  0 

•13 

-10 

.01 

-20  j 

5.  From  figure  8,  it  is  determined  that  the  factor  by  which  to 
multiply  the  liquid-water-content  values  presented  in  step  4  is  1.3. 

This  particular  value  corresponds  to  the  design  value  of  cloud  horizontal 
extent  suggested  in  figure  8  for  class  I  (instantaneous)  icing  conditions 
with  Pe  taken  as  0.001.  A  tabulation  of  the  final  liquid-water -content 
values  suitable  for  design  purposes  is  presented  in  the  following  table: 

Instantaneous  Conditions  Corresponding  to 
Pe  =  0.001 

[Horizontal  extent  considered] 


w 

T 

D 

3.44 

32 

3.11 

10 

2.48 

0 

15 

1.56 

-10 

.71 

-20 

1.95 

32 

1.32 

10 

1.17 

0 

30 

.63 

-10 

.23 

-20 

.74 

32 

•63 

10 

.41 

0 

45 

.17 

-10 

.01 

-20 

2k 
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In  the  above  table,  care  should  be  taken  to  note  that  the  exceedance 
probability,  Pe  =  0.001,  applies  to  each  icing  condition  listed  and  not 
to  the  class  in  general. 


The  Partial  Probability  Charts 


The  partial  probability  charts  essentially  provide  a  rapid  means  of 
determining  the  approximate  probability  of  exceeding  any  specified  group 
of  values  of  liquid-water  content  associated  simultaneously  with  tempera¬ 
ture  and  drop-diameter  values  lying  within  specified  ranges.  Thus,  if 
an  anti -icing  component  (e.g.,  unit  span  of  a  jet -inlet  guide  vane)  is 
designed  by  the  method  of  reference  16  to  withstand  a  constant  value6 
of  weight  rate  of  water-drop  impingement  Mt,  as  computed  from  a  partic¬ 
ular  set  of  values  of  liquid-water  content,  temperature,  and  drop  diam¬ 
eter,  and  if  other  combinations  of  these  same  variables  are  found  which 
correspond  to  the  same  weight  rate  of  drop  impingement,  then  the  prob¬ 
ability  of  encountering  conditions  under  which  the  component  will  not 
perform  satisfactorily  can  be  determined  approximately.  The  calculation 
of  this  probability  can  be  performed  with  the  aid  of  the  partial  prob¬ 
ability  charts  (fig.  10). 


General  procedure  for  using  the  partial  probability  charts.-  To  use 
the  partial  probability  charts,  the  following  procedure  is  suggested: 


1.  Establish  values  of  liquid-water  content,  Wj_,  for  various  com¬ 
binations  of  drop  diameter  and  temperature,  which  are  regarded  as  crit¬ 
ical  (a  marginal  surface)  for  a  given  component  of  a  thermal  system. 

Values  of  Wj_  can  be  obtained  by  computing  values  of  liquid-water  content 
required  to  provide  a  constant  weight  rate  of  water  impingement  on  the 
component .  The  equation 

w  =  - — -  (5) 

1  0.33EmVtc 


®The  choice  of  a  constant  value  of  weight  rate  of  water-drop  impingement 
rather  than,  for  example,  the  use  of  the  amount  of  heat  to  raise  the 
surface  temperature  to  32°  E,  provides,  according  to  reference  l6,  a 
reasonable  basis  for  estimating  the  probability  of  encountering  an 
icing  condition  too  severe  for  the  ice -protection  equipment.  The 
amount  of  heat  supplied  for  the  particular  weight  rate  of  water  impinge¬ 
ment  selected  depends  upon  the  type  of  protection  required  for  the 
component,  that  is,  whether,  for  example,  all  the  water  impinging 
should  be  evaporated  instead  of  having  the  surface  temperature  merely 
brought  to  32°  F.  Since  for  various  water-drop  sizes  and  temperatures 
the  amount  Of  heat  required  for  protection  is  not  unique  for  a  given 
weight  rate  of  water-drop  impingement,  the  marginal  surface  will  not 
represent  the  same  degree  of  protection  for  all  combinations  of  drop 
size  and  temperature.  However,  it  will  very  nearly  do  so.  As  a  result, 
the  probability  value  computed  from  use  of  the  marginal  surface  is  not 
exactly  correct. _ _ _ 


derived  in  reference  l6 ,  can  be  used  to  calculate  the  values  of  as 

a  function  of  drop  diameter  and  temperature,  since  Em  is  a  function 
of  these  latter  two  variables.  The  particular  value  of  Mt  used  to 
calculate  values  of  Wi  is  obtained  by  inserting  into  equation  (5) 
initial  values  of  Wi,  Em,  V,  t,  and  c  corresponding  to  a  selected 
meteorological  design  condition  and  airfoil  section.  The  appropriate 
values  of  Ej-^  to  employ  in  equation  (5)  for  the  chosen  values 
of  V,  t,  and  c  must  be  determined  by  using  some  type  of  water-drop- 
trajectory  calculation,  such  as  presented  in  reference  13 .  When  using 
equation  (5)  to  obtain  values  of  W±,  the  computations  should  be  made 
for  all  the  possible  combinations  of  Tj_  and  D^  listed  in  table  II. 

In  total,  there  are  such  combinations. 

2.  Determine  the  probability,  APj_,  of  exceeding  each  of  the  values 
of  liquid-water  content,  Wi,  determined  in  step  1.  To  determine  a 
specific  value  of  APi  from  a  particular  partial  probability  chart 
predented  in  figure  10,  the  procedure  is  as  follows:  Follow  a  vertical 
line  representing  the  liquid-water  content  until  the  line  intersects  a 
horizontal -extent  line.  From  the  horizontal-extent  line,  follow  hori¬ 
zontally  until  an  intersection  with  a  temperature-interval  contour  is 
obtained.  From  this  Intersection,  follow  a  vertical  ~Hnp  downward  to 
obtain  an  intersection  with  a  diameter-interval  line.  A  horizontal  line 
through  this  latter  intersection  yields  the  value  of  AP*  on  the  ordi¬ 
nate  scale.  When  it  happens  that  a  horizontal  line  drawn  from  the 
horizontal-extent  curve  will  not  intersect  the  desired  temperature- 
interval  line,  all  that  can  be  said  about  the  value  of  APi  for  a  given 
diameter  interval  Is  that  it  is  smaller  than  the  value  the  horizontal 
line  would  indicate  when  projected  to  the  extreme  left  side  of  the  par¬ 
tial  probability  chart.  Such  values  are  usually  Bmall  enough  to  be 
neglected. 

3.  Establish  the  probability,  P,  of  exceeding  all  the  values  of 
liquid-water  content  defined  in  step  1.  This  probability  may  be  obtained 
by  summing  all  the  individual  values  of  APi,  determined  in  step  2,  in 
accordance  with  equation  (3).  The  number  of  APi  values  which  must  be 
summed  is  k2,  since  this  number  represents  all  combinations  of  tempera¬ 
ture  and  drop  diameter  possible  by  using  the  charts  of  figure  10. 

Example  showing  the  use  of  the  partial  probability  Charts.-  The 
example  of  the  jet-engine  inlet-guide  vane  will  be  continued  to  demon¬ 
strate  the  use  of  the  partial  probability  charts  for  a  vane  section 
assumed  to  be  heated  to  100-percent  chord.  The  starting  point  of  this 
example  will  be  taken  with  established  sets  of  values  of  liquid-water 
content,  drop  diameter,  and  temperature  that  define  a  marginal  surface 
for  a  unit  span  of  the  vane.  These  sets  of  values,  which  are  presented 
in  table  VI,  will  be  presumed  to  have  been  obtained  by  calculating 
values  of  W±  from  equation  (5),  using  a  particular  design  icing  condi¬ 
tion  from  the  equiprobability  charts.  The  particular  icing  condition 
presumed  to  have  been  selected  for  the  computations  is  one  taken  from 
the  example  use  of  the  equiprobability  charts  and  defined  by  the  following 
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values:  ¥  =  2.48  grams  per  cubic  meter;  T  =  0°  F;  D  =  15  microns;  and 
cloud  horizontal  extent,  0.5  mile.  With  these  available  data,  the  steps 
employed  to  use  the  partial  probability  charts  are: 

1.  Values  of  APi  are  selected  from  the  marginal  probability 
Charts  for  Pacific  coast  cumulus  (fig.  10(b))  to  correspond  with  the 
conditions  chosen  from  the  equiprob ability  chants  for  a  horizontal 
extent  of  0.5  mile.  The  results  are  tabulated  in  table  VI  for  the  4-2 
different  possible  combinations  of  Tj_  and  D-j_. 

2.  The  probability  of  encountering  an  icing  condition  lying  out¬ 
side  the  surface  defined  in  W-T’-D  space  by  the  Wi-Tj_-Dj_  values 
listed  in  table  VI  is  obtained  by  adding  the  partial  probabilities  found 
for  the  various  conditions  presented  in  the  table.  In  the  summation, 
values  entered  in  table  VI  as  being  less  than  a  certain  specified  amount 
are  neglected  because  the  values  are  exceedingly  small  and  not  exactly 
determined  in  magnitude.  These  values  arise  from  the  fact  that  the 
charts  are  not  constructed  to  encompass  quite  small  values  of  AP^.  For 
the  example  presented  in  table  VI,  the  value  of  the  probability,  P,  was 
found  to  be  0 . 0657 •  The  value  of  0.0657  can  be  interpreted  to  mean  that 
in  approximately  93  out  of  100  icing  encounters  the  values  of  liquid- 
water  content,  for  all  combinations  of  drop  diameter,  will  lie  below  the 
marginal  surface  and  hence  be  encompassed  by  the  design. 


CONCLUDING  REMARKS 


The  equiprobability  surfaces  and  marginal  probability  calculation 
charts  presented  herein  provide,  for  the  United  States,  a  representation 
of  the  data  now  available  on  the  meteorological  factors  responsible  for 
aircraft  icing,  expressed  in  terms  of  the  probability  of  exceeding 
various  combinations  of  values  of  these  factors.  The  usefulness  of  the 
results  is  limited  somewhat  by  the  amount  of  data  available;  errors  of 
measurement,  especially  at  large  values  of  drop  diameter;  nonrepresent¬ 
ativeness  of  the  flight  procedures  during  research  flights;  and  limited 
altitude  and  geographical  extent  of  the  research  flight  program. 

In  spite  'of  these  limitations,  the  probability  analysis  presented 
does  provide  an  indication,  heretofore  unavailable,  of  the  combinations 
of  icing  conditions  having  equal  probability  of  being  exceeded  in  the 
United  States,  and  also  the  probability  of  exceeding  a  special  set  of 
icing  conditions.  Also,  a  procedure  has  been  established  for  the  sta¬ 
tistical  analysis  of  future  icing  meteorological  data  obtained  on  a 
world-wide  basis. 


Ames  Aeronautical  Laboratory 

National  Advisory  Committee  for  Aeronautics 

Moffett  Field,  California,  April  18,  1952 
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APPENDIX  A 

DETAILS  OF  DETERMINING  THE  EQUIPROBABILITY 
SURFACE  CONTOURS 


The  distribution  functions  used  in  the  equiprobability  study  are 
the  normal  and  the  Gurnbel  probability  distributions.  Graphical  methods 
were  used  to  evaluate  the  distributions  by  using  normal -probability 
graph  paper  for  the  normal  distributions  and  a  specially  constructed 
paper  for  the  Gumbel  distributions.  These  special  types  of  graph  paper 
have  the  property  that  the  cumulative  distribution  curves  are  repre¬ 
sented  by  straight  lines. 

The  Gumbel  distribution,  taken  from  reference  11,  is  given  by  the 
equation 

PG  =  l-ee  (Al) 

In  equation  (Al),  the  constant  u  is  the  mode  and  a  is  a  quantity 
which  measures  the  concentration  of  the  frequency  distribution  about  the 
mode.  Since  on  the  specially  constructed  Gumbel  distribution  graph  paper 
equation  (Al)  will  appear  as  a  straight  line,  it  is  desirable  to  define 
the  line  by  specifying  two  points  on  the  curve.  For  the  purposes  of 
this  analysis,  the  points  are  chosen  at  which  the  probability  Pq  has 
the  values  of  O.63  and  0.05*  The  value  of  Pq  =  O.63  may  be  found  by 
placing  x  =  u  and  the  value  of  Pq  =  0.05  may  be  found  by 
placing  a(x-u)  =3.0. 

When  a  straight-line  fit  could  not  be  obtained  from  the  data  by 
using  Gumbel  distribution  paper,  an  arbitrary  smooth  curve  on  normal 
probability  paper  was  used  to  represent  the  distribution  function.  This 
procedure  conserved  the  essential  form  of  the  observed  distribution, 
smoothed  out  the  irregularities,  and  provided  a  reasonable  basis  for 
limited  extrapolation. 

Cases  2  and  5  (table  I)  are  used  to  demonstrate  how  the  normal  and 
Gumbel  probability  distributions  were  applied  to  the  data  since  these 
cases  were  the  ones  in  which  the  largest  number  of  encounters  with  icing 
were  recorded.  As  mentioned  in  the  text,  the  general  procedure  is  to 
evaluate  the  functions  Q^),  R(W,T'),  S(D,W,T>),  R'(W,D,T') 
and  S»(D,T‘). 

Case  2,  Cumulus  Clouds  in  the  Pacific  Coast  Area 

The  distribution  functions.-  The  distribution  of  the  temperature 
depression  below  freezing  of  the  most  severe  icing  condition  in  each 
icing  encounter  for  this  case  is  shown  in  figure  11  plotted  on  Gumbel' s 
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distribution  paper.  In  this  figure  the  broken  line  represents  the 
observed  values  and  the  straight  line  is  the  estimated  line  of  best  fit. 
This  line  defines  the  function  Q(T’). 

The  distribution  of  maximum  liquid-water  content  per  encounter  was 
plotted  for  all  values  of  T'  equal  to  or  greater  than  each  of  the 
following  values:  0,  8,  12,  17,  22,  27,  28,  32,  and  37-  The  plots 
for  T'>0,  and  T’>27  are  presented  in  figure  12.  The  values  of  the 
distribution  parameters,  W0 i63  and  W0>05,  obtained  from  these  distribu¬ 
tions  and  plotted  as  functions  of  T',  are  shown  in  figure  13.  The 
smoothed  curves  in  figure  13  define  the  liquid-water-content  distribu¬ 
tion  function,  R(W,T').  These  curves  were  drawn  to  the  spontaneous 
freezing  temperature,  T'  =  72  F°. 

The  distribution,  S(D,W,T')  of  drop  size  to  be  used  with  the  dis¬ 
tribution  of  maximum  liquid-water  content  per  encounter  is  to  be  deter¬ 
mined  for  all  observations  with  W  and  T'  greater  than  certain  values. 
Small  values  of  W  need  not  be  considered  in  this  case  since  conditions 
of  small  W  and  large  D  are  to  be  expressed  in  terms  of  S’  and  R'. 

The  distribution,  S(D,W,T’)j  of  drop  diameter  for  all  values  of  T’ >0 
and  W>0.5  is  shown  in  figure  l4  plotted  on  normal  probability  paper. 
Similar  plots  were  made  for  W>0.5,  T'  >22,  and  WlkO.8,  T'^.0;  and 
since  these  distributions  were  very  similar  to  the  one  for  W^.0-5 , 

T' >0,  a  statistical  test,  the  chi-square  test  (reference  17)  was 
applied  to  test  the  significance  of  the  differences.  The  results  of 
this  test  showed  that  the  differences  between  these  distributions  were 
no  greater  than  would  be  expected  as  a  result  of  random  sampling;  hence 
it  was  concluded  that  any  dependence  of  S  upon  ¥  and  T'  for  values 
of  W>0.5  was  too  small  to  be  reliably  indicated  by  the  data.  Accord¬ 
ingly.,  the  distribution  S(D,W,T' )  shown  in  figure  l^was  used  for  all 
values  of  W>0.5. 

In  order  to  define  the  equiprobability  surfaces  in  the  region 
where  W  is  small  and  D  is  large,  the  distributions  S'(D,T') 
and  R'(W,D,T')  must  be  determined.  The  distribution,  ST(D,T’)>  of 
the  maximum  value  of  D  per  icing  encounter  was  examined  for  T'^.0, 

T*  ^22,  and  T’>32;  and  the  application  of  the  chi-square  test  to  these 
distributions  indicated  that  there  was  no  significant  detectable  depend¬ 
ence  of  S *  on  temperature.  The  distribution,  S'(D,T')j  for  T’>0  is 
shown  in  figure  15- 

In  order  to  determine  R’(W,D,T' ),  the  distribution  of  W  for  large 
values  of  D,  this  distribution  was  plotted  for  the  following  cases: 

D>20,  T'>0;  D  —  30,  T' >0;  D>44,  Tf>0;  D>20,  T’>22;  and  D >  30,  T’>22. 
A  study  of  these  distributions  indicated  that  there  was  apparently  no 
effect  of  T'  on  the  variation  of  R'  with  W.  There  was,  however, 
a  significant  effect  of  D  on  the  variation  of  R'  with  W. 

Figure  16  shows  the  distributions  for  T*  >0  and  D^20,  30,  and  44. 

Figure  17  shows  the  parameters  defining  the  R'  distributions  for 
T*  >0  plotted  as  a  function  of  D.  It  is  likely  that  the  failure 
of  the  data  to  show  for  large  values  of  D  a  decrease  of  W  with 
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increasing  T»,  as  was  the  case  in  figure  13,  is  due  to  the  fact  that 
insufficient  data  are  available  for  large  values  of  D  to  define  the 
function  R*  when  TT  is  greater  than  22,  since  it  is  only  in  this 
range  that  the  function  R  becomes  strongly  dependent  upon  T*.  In 
the  absence  of  more  complete  information,  it  has  been  assumed  that 
when  T'  ->22,  the  dependence  of  R'  on  T'  is  similar  to  the  dependence 
of  R  on  T' .  Accordingly,  the  values  of  the  R'  distribution  param¬ 
eters  (the  mode  W0#63  and  the  fifth  percentile  WQr05)  shown  in  fig¬ 
ure  17  are  altered  to  obtain  values  of  the  distribution  parameters  which 
are  considered  applicable  when  T'  >22.  x  The  factor  used  to  alter  any 
particular  value  of  W  presented  in  figure  17  is  established  from  fig¬ 
ure  13  by  forming  the  ratio  of  the  value  of  W  at  a  selected  value 
of  T*  to  the  value  of  W  applicable  to  values  of  T' <22.  The  values 
of  .63  wo  .05  which  define  the  function  R'  were  obtained  in 
this  manner,  and  are  presented  in  table  VII. 

Construction  of  the  equiprobability  contours.-  After  the  distribu¬ 
tion  functions  have  been  established,  the  equiprobability  contours  for 
various  temperature  depressions  can  be  determined.  The  following 
examples  illustrate  the  use  of  the  probabilities  Q,  R,  S,  R',  and  S’ 
in  determining  the  equiprobability  surface  for  Pacific  coast  cumulus 
clouds  corresponding  to  an  exceedance  probability  of  0.01.  For  a  value 
of  Pe  =  0.01,  the  equation 

'  QRS  =  0.01  (A2) 

applies  to  large  values  of  W.  Similarly,  the  equation 

QR'S*  =  0.01  (A3) 

applies  to  small  values  of  W.  In  evaluating  equations  (A2)  and  (A3), 
two  equiprobability  surface  contours  are  selected.  These  contours  are 
for  T’  =  0  and  T'  =32. 

Large  values  of  W:  Table  VIII  shows  the  calculation  of  W, 
using  equation  (A2)  for  various  arbitrarily  chosen  values  of  D.  Com¬ 
putations  are  shown  when  T’  =  0  and  when  T*  =32. 

When  T’  =  0,  it  should  be  noted  that  Q  =  1.0.  Thus,  equa¬ 
tion  (A2)  becomes 

R  =  0.01/S  (A4) 

The  values  of  S  necessary  in  equation  (Ah)  to  calculate  values  of  R 
for  various  values  of  D  are  obtained  from  figure  l4.  Values  of  R 
calculated  from  equation  (A4)  are  used  to  define  the  desired  values 
of  W  on  the  curve  depicting  the  distribution  function  R(W,T’ )  for 
a  value  of  T*  =  0.  This  function  is  defined  as  a  straight  line  on 
Gumbel  distribution  paper  by  establishing  the  values  of  W  for  the 
mode  and  at  the  fifth  percentile.  The  values,  W0^63and  W0  #05,  used 
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in  defining  the  straight-line  distribution  of  R(W,T*)  are  obtained 
from  figure  13  corresponding  to  a  value  of  T*  =  0.  For  this  partic¬ 
ular  case,  the  result  would  be  the  same  as  presented  in  figure  12. 

The  desired  values  of  W  chosen  from  the  straight-line  distribution 
for  the  calculated  values  of  R  are  presented  in  table  VIII,  part  (a). 

When  T'  =  32,  Q  is  found  from  figure  11  to  have  the  value 
of  0.176.  Hence,  the  equation  expressing  values  of  R  is 


0.01 

0.176S 


(A5) 


Values  of  S  required  to  solve  equation  (A5)  are  obtained  again  from 
figure  lk  for  various  arbitrarily  selected  drop  diameters.  To  find 
values  of  W  to  which  the  calculated  values  of  R  correspond,  it  is 
again  necessary  to  establish  R(W,T')  as  a  straight-line  function 
of  W  on  Gumbel  distribution  paper.  The  two  points  required  to 
establish  the  straight-line  variation  are  the  values  of  W  at  the 
mode  (W0#63)  and  at  the  fifth  percentile  (W0>05).  These  values,  for 
a  value  of  T'  =  32,  are  found  from  figure  13  to  be  WQ >63  =  0.28 
and  W0#05  =  1.26.  The  values  of  W  found  from  the  straight-line 
variation  of  W  with  R(D,T* )  are  tabulated  in  table  VIII,  part  (b), 
along  with  the  corresponding  values  of  D,  S,  and  R. 

Small  values  of  W:  Sections  of  the  equiprobability  contours 
corresponding  to  small  values  of  W  are  calculated  in  identically  the 
same  manner  as  for  large  values  of  W.  In  both  the  case  of  T*  =  0 
and  of  T'  =  32,  values  of  S'  were  obtained  from  figure  15  for  cal¬ 
culating  values  of  R*  from  the  equation 

R'  =  P/QS? 

The  values  of  W  at  the  mode  and  the  fifth  percentile  for  defining  W 
as  a  function  of  R*(W,D,T')  were  obtained  from  table  VII.  The  calcu¬ 
lations  for  T*  =  0  are  presented  in  table  IX,  part  (a),  and  the  cal¬ 
culations  for  T*  =  32  are  presented  in  table  IX,  part  (b). 

Intermediate  values  of  W:  From  tables  VIII  and  IX,  equiprobabil¬ 
ity  curves  can  be  plotted  for  both  high  and  low  values  of  W,  as  is 
shown  in  figure  18.  For  the  intermediate  values  of  W,  however,  the 
final  equiprobability  contours  must  be  obtained  by  joining  the  upper 
and  the  lower  portions  by  a  smooth  connecting  line,  shown  by  a  broken 
line  in  figure  18.  The  upper  portions  of  the  final  contours  represent 
maximum  values  of  liquid-water  content  per  icing  encounter  and  the 
corresponding  drop  diameter;  the  lower  portions  represent  maximum  drop 
size  per  encounter  and  the  corresponding  liquid-water  content;  and  the 
intermediate  portions  may  represent  either  one  of  these  combinations. 
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The  equiprobability  contours  presented  in  figure  18  are  an 
incomplete  set  of  those  presented  in  figure  4(b). 


Case  5,  Layer  Clouds  in  Eastern  United  States 


The  distribution  functions.-  The  distribution  of  the  temperature  of 
the  most  severe  icing  observation  in  each  icing  encounter  for  this  case 
is  shown  in  figure  19.  Since  a  satisfactory  straight-line  representa¬ 
tion  could  not  be  found  for  this  frequency  distribution,  the  smooth 
curve  on  normal  probability  paper  given  in'  figure  19  waB  used  to  define 
the  function  Q(TT)« 

The  distribution  of  maximum  liquid-water  content  per  icing 
encounter,  R(W,T’),  was  plotted  for  values  of  Tf  equal  to  or  greater 
than  0,  7,  12,  17,  22,  and  27.  The  curves  for  Tf  >0  and  T*  «£.17  are 
presented  in  figure  20.  The  values  of  WQ-63  and  WD  Q5  obtained  from 
each  of  the  six  distribution  curves  were  plotted  as  functions  of  T* 

(fig.  21),  and  smooth  curves  were  drawn  to  represent  the  distribution 
parameters  which  define  the  function  R(w,T* ) . 

The  distribution  of  drop  diameter  for  values  of  W>0.2  and  T’  >0 
is  shown  in  figure  22.  As  in  case  2,  this  distribution  was  found  to  be 
approximately  independent  of  W  and  T'  for  values  of  ¥  greater 
than  0.2.  The  curve  in  figure  22  was,  therefore,  used  to  represent  the 
function  S(D,W,T’). 

The  distribution  of  the  largest  value  of  D,  S’(D,T)  observed  in 
each  icing  encounter  is  shown  in  figure  23.  This  distribution  curve, 
which  was  found  to  be  approximately  independent  of  T* ,  was  used  to 
define  the  function  S’(D,T’). 

In  order  to  determine  R'(¥,D,T')>  the  distribution  of  W  corre¬ 
sponding  to  large  values  of  D,  the  distribution  of  W  was  plotted  for 
all  values  of  T*  and  for  values  of  D  equal  to  or  greater 
than  15,  17,  19>  21,  23,  and  25-  Values  of  ¥0>63  and  WQ  >05  for  each 
of  these  distributions  were  plotted  against  D  ‘as  Shown  in  figure  24. 

As  a  result  of  the  small  amount  of  data  available  for  large  values 
of  D,  there  is  a  considerable  scatter  in  these  distribution  parameters. 
The  smooth  curves  shown  in  figure  24  were  drawn  to  represent  approxi¬ 
mately  the  dependence  of  R '  upon  D.  To  determine  the  dependence 
of  R*  upon  T’,  cases  of  D>17  were  chosen  since  this  was  the 
largest  value  of  D  for  which  sufficient  data  were  available  to  provide 
a  reasonably  reliable  sample.  Distributions  of  W  for  D>  17  and  T’>0, 
12,  and  17  were  plotted  and  the  resulting  values  of  W0.63  and  W0.os 
from  these  distributions  were  used  with  the  data  from  figure  24  to  pro¬ 
vide  the  basis  for  the  construction  of  a  set  of  curves  defining  the 
parameters  of  the  distribution,  R»(W,D,T»),  (fig.  25).  Curves  repre¬ 
senting  values  of  WQ  ^3  and  WQ  >Q5  for  D  >17  were  drawn  through  the 


32 


NACA  TN  2738 


data  points  (fig.  25)#  using  the  shape  of  the  curves  of  figure  21  (the 
distribution  parameters  for  R)  as  a  guide.  The  curves  for  other  values 
of  D  were  then  drawn  from  the  points  determined  by  the  values 
for  T*  =0  (obtained  from  fig.  24)  in  such  a  way  as  to  conform  to  the 
pattern  established  by  the  curves  for  D>17*  These  curves  define  the 
function  R*(W,D,T'). 

Values  of  the  probability  terms,  Q,  R,  S,  R',  and  S',  were  obtained 
with  the  aid  of  figures  19,  20,  22,  25,  and  23,  respectively,  and  were 
used  to  construct  the  equiprob ability  surfaces  for  values  of  Pe 
from  0.1  to  0.001,  which  are  presented  in  figure  7* 


Cases  1,  3,  and  4 


By  following  procedures  similar  to  those  described  in  the  preceding 
examples,  equiprobability  surfaces  also  could  be  determined  for  the 
following  cases:  (1)  layer  clouds  in  the  Pacific  coaBt  area  (fig.  3)j 
(2)  layer  clouds  in  the  plateau  area  (fig.  5)j  and  (3)  cumulus  clouds  in 
the  plateau  area  (fig.  6).  The  case  of  cumulus  clouds  in  the  eastern 
United  States  was  omitted  because  of  insufficient  data. 
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APPENDIX  B 

DETAILS  OF  THE  PREPARATION  OF  THE  PARTIAL 
PROBABILITY  CHARTS 


The  same  eases  as  were  used  in  describing  the  construction  of  the 
equiprobability  charts  will  be  used  to  describe  in  detail  the  construc¬ 
tion  of  the  partial  probability  charts.  The  general  procedure  in  con¬ 
structing  the  charts  is  first  to  evaluate  the  equation 

APi  *  PAT  PAD 

The  resulting  values  of  AP^  for  different  combinations  of  liquid-water 
content  and  intervals  of  drop  diameter  and  temperature  are  then  arranged 
into  a  graphical  form  for  convenient  use. 


Evaluation  of  AP-^ 


Case  5,  layer  clouds  in  eastern  United  States.-  Values  of  PAt  for 
this  case  were  obtained  from  figure  19  which  provides  information  as  to 
the  probability  that  the  temperature  depression  will  be  lower  than  any 
particular  value.  To  find  a  value  of  Pat  (the  probability  that  the 
temperature  will  lie  within  the  interval  AT)  from  figure  19,  it  is  only 
necessary  to  take  the  difference  between  the  probabilities  indicated  for 
the  two  chosen  limits  of  AT.  Values  of  PAt  determined  from  figure  19 
are  presented  in  table  X,  part  (a),  for  selected  temperature  intervals. 

Values  of  Pad  (the  probability  that  a  drop  diameter  will  lie  within 
the  interval  AD)  were  obtained  in  a  manner  similar  to  that  used  in 
obtaining  values  of  PAt«  In  this  case,  however,  a  composite  drop- 
diameter  distribution  curve  was  used  in  order  to  allow  approximately  for 
the  occurrence  of  icing  encounters  in  which  the  conditions  of  maximum 
severity  are  determined  primarily  by  drop  size  rather  than  liquid-water 
content.  The  composite  distribution  curve  is  presented  in  figure  2 6 
where  the  probability  values  in  this  special,  case  are  called  S'*.  As 
can  be  seen  from  figure  26,  the  composite  distribution  curve  follows 
the  distribution  of  all  observed  values  of  drop  diameter  for  the  lowest 
90  percent  of  the  distribution  and  the  distribution  of  max-1  mum  observed 
value  per  encounter  (fig.  23)  for  the  highest  1  percent  of  the  cases. 

The  intermediate  9  percent  of  the  cases  is  represented  by  a  straight 
line  connecting  the  two  distributions.  Values  of  Fad  determined 
through  use  of  figure  26  are  presented  in  table  X,  part  (b),  for 
selected  drop-diameter  intervals. 

In  order  to  determine  values  of  Ify (AT),  first  the  distributions 
of  maximum  liquid-water  content  per  encounter  were  plotted  on  Gurnbel 
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distribution  graph  paper  for  each  of  the  chosen  temperature  intervals 
above  -9* 5°  F  in  a  manner  similar  to  that  done  in  figure  20  for  two 
ranges  of  values  of  T'.  Values  of  W0  ,63  and  W0  >05  thus  could  be 
obtained;  and  these  values  were  plotted  as  a  function  of  temperature 
and  extrapolated  to  obtain  values  of  W0  ,63  and  WD ,os  applicable  to 
temperature  intervals  below  -9.5°  F.  The  resulting  values  of  WQ>63 
and  W0.05  for  all  temperature  intervals  are  presented  in  table  X, 
part  (c).  In  addition,  there  are  presented  values  of  W  which  were 
chosen  arbitrarily  from  the  straight-line  distribution  established  by 
the  two  points,  W0.63  and  W0>05.  Corresponding  to  the  selected  values 
of  W  exist  values  of  Pv^(AT)  which  also  can  be  read  from  the  straight- 
line  distribution  defined  by  the  points  W0  >63  and  W0  <05 .  These  values 
of  Py^AT)  are  also  presented  in  table  X,  part  (c). 

After  the  values  of  P^p,  P^p,  and  Py^  (AT )  are  determined,  values 
of  APi  can  be  determined  simply  by  multiplying  the  constituent  prob¬ 
abilities  together  for  various  combinations  of  liquid-water  content  and 
intervals  of  temperature  and  drop  diameter.  Such  computations  are  sum¬ 
marized  in  table  XI  for  two  temperature  intervals  (32*0°  to  20.5°  F 
and  20. 5°  to  10.5°  F),  two  drop-size  intervals  (0  to  9*5  microns 
and  9*5  to  12. 5  microns)  and  a  range  in  liquid-water-content  values. 

The  values  of  APi  listed  in  the  table  are  ones  which  can  be  obtained 
through  use  of  the  partial  probability  charts. 

Other  caseB.-  The  quantities  used  in  calculating  values  of  APi 
for  the  remaining  cases  studied  were  obtained  in  the  same  manner  as 
case  5,  with  one  exception.  In  case  2,  the  distribution  of  the  greatest 
drop  diameter  per  encounter  was  used  for  the  upper  5  percent  of  the 
diameter  distribution  instead  of  the  upper  1  percent.  The  reason  for 
this  deviation  was  because  of  the  greater  frequency  Of  encounters  in 
which  the  most  severe  icing  condition  waB  determined  by  the  maximum 
value  of  drop  size. 


Example  Construction  of  the  Partial  Probability  Charts 


Rather  than  compute  values  of  APj_  every  time  they  are  required, 
use  of  the  data  presented  in  tables  X  and  XI  in  chart  form  is  a  great 
convenience.  In  general,  the  technique  of  constructing  a  suitable  chart 
for  evaluating  values  of  APi  utilizes  the  relation 

log  APi  =  log  P^p  +  log  +  log  Ifyi(AT)  (Bl) 

in  such  a  manner  that  the  chart  performs  the  function  of  adding  the 
right-hand  terms  in  the  equation. 

The  various  steps  in  the  construction  of  a  chart  to  solve  equa¬ 
tion  (Bl)  (derived  from  equation  (4))  are  shown  in  figure  27*  Fig¬ 
ure  27(a)  establishes  the  liquid-water  content  as  a  function  of  the  sum 
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of  the  logarithms  of  P&r  and  Pw±(AT).  Since  the  value  of  P^  corre¬ 
sponding  to  a  given  temperature  interval  is  a  constant,  the  liquid-water- 
content  distribution  curve  is  displaced  horizontally  from  the  origin  by 
the  amount  of  log  P^.  The  values  of  P^,  PWi(Zfl’)  and  used  in 
preparing  figure  27(a)  were  obtained  from  table  X,  parts  (a)  and  (c), 
for  a  temperature  interval  of  32°  to  20.5°  F. 

The  next  step  in  the  construction  is  shown  in  figure  27(b)  where 
the  liquid-water-content  scale  on  the  right  is  transformed  into  a  corre¬ 
sponding  scale  at  the  top  by  use  of  a  line  drawn  at  an  angle  of  45°  on 
the  grid  system.  This  line  may  be  labeled  the  10-mile  horizontal-extent 

because  the  ratio  of  any  value  on  the  top  scale  to  the  correspond¬ 
ing  value  on  the  right-side  scale  is  unity.  For  example,  in  the  figure 
a  value  of  Wj[  =  0.5  on  the  top  scale  leads  to  a  value  of  w^  =  0.5 
on  the  vertical  scale  when  the  10-mile  horizontal-extent  n™>  is  used. 

For  other  horizontal  extents,  however,  a  factor  must  be  applied  to  the 
liquid-water  content  for  10-mile  horizontal  extent  to  obtain  the  value 
appropriate  to  the  extent  under  consideration.  Accordingly,  another 
line  (derived  from  reference  3)  has  been  drawn  on  figure  27(b)  to  incor¬ 
porate  automatically  the  factor  which  should  be  applied  to  a  50-mile 
horizontal  extent  at  different  values  of  PAT  and  PWi(/d1).  Thus,  a 
value  of  =  0.5  on  the  top  scale  would  be  transformed  into  a  value 
of  Wi  =  0.82  on  the  right-hand  scale  for  a  horizontal  extent  of 
50  miles.  The  point  bn  the  liquid-water  distribution  curve  to  which  the 
value  of  W±  =  0.82  applies  is  labeled  point  A  in  figure  27(b). 

The  only  other  factor  which  must  be  included  in  the  partial  prob- 
charts  is  the  probability  factor  P/vn»  A  method  of  incorporat¬ 
ing  the  term  is  by  using  a  scale-transformation  reference  line  (similar 
in  operation  to  the  horizontal-extent  lines).  Figure  27(c),  the  top 
part  of  which  is  the  same  as  figure  27(b),  is  presented  to  show  the 
function  of  the  scale-transformation  reference  line  in  determining  the 
total  sum  of  log  Pat*  log  (AT ) ,  and  log  PAt).  As  an  example,  a  value 
of  Wi  =  0.5  is  taken  for  which  it  is  desired  to  find  the  sum  of  the 
logs  of  the  three  probabilities,  P^j,  PWl(M'),  and  P^p  applicable  to 
a  10-mile  horizontal  extent,  a  temperature  interval  of  32°  to  20.5°  F 
and  a  drop-diameter  interval  of  0  to  9.5  microns.  Accordingly,  a  value 
of  Wi  =  0-5  is  selected  on  the  liquid -water-content  scale  and  is  fol¬ 
lowed  downward,  vertically,  until  an  intersection  with  the  10-mile 
horizontal-extent  line  is  obtained  (point  1).  From  point  1,  a  line  is 
followed  horizontally  until  an  intersection  with  the  liquid-water-content 
distribution  curve  for  the  temperature  interval  of  32°  to  20.5°  F 
(point  2).  The  horizontal  distance  that  point  2  is  away  from  the  log 
reference  line  is  a  measure  of  the  sum  of  log  PArp  and  log  P^.  (AT).  If 
a  vertical  line  from  point  2  is  followed  downward  until  an  intersection 
with  the  scale  transformation  reference  line  is  reached  (point  3),  the 
absolute  value  of  the  sum  of  log  Pat1  and  log  Pj^ (AT)  can  be  read  on  the 
Beales  either  to  the  left  or  to  the  right  of  the  point.  However,  if  the 
vertical  line  followed  from  point  2  is  continued  downward  until  the 
intersection  with  the  drop-diameter  interval  line  is  reached  (point  4), 
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the  sum  of  the  logs  of  Py^AT)  and  can  be  read  from  the  left- 

hand  scale  for  log  of  APi.  The  reason  why  the  total  can  be  read  is 
that  the  drop-diameter  interval  line  is  displaced  vertically  downward 
from  the  scale  transformation  reference  line  by  the  amount  of  log  P^p 
for  the  diameter  interval  of  0  to  9*5  microns.  The  actual  value 
of  Pap  used  in  the  construction  was  obtained  from  table  X,  part  (b ) . 

The  final  value  of  log  APi  read  from  the  partial  probability 
chart  in  figure  27(c)  is  the  result  sought,  namely,  a  solution  of  equa¬ 
tion  (Bl).  For  most  rapid  use  of  the  chart,  however,  it  is  desirable 
to  have  values  Of  APi  obtainable  directly.  This  added  convenience 
can  be  incorporated  in  a  chart  by  making  each  integer  of  the  logarithmic 
scales  in  figure  27(c)  correspond  to  a  cycle  on  semilogarithmic  graph 
paper.  The  partial  probability  charts  presented  in  figure  10  have  been 
constructed  to  include  this  feature,  and  hence  will  yield  values 
of  APi  directly. 
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TABLE  I.-  CLASSIFICATION  OF  ICING  CONDITIONS 
AND  AMOUNT  OF  DATA  AVAILABLE  FOR  EACH  CASE 


Case 

Geographical 

area 

Cloud 

type 

Number  of  icing 
encounters 

Number  of  observations 
of  liquid-water  content 
and  drop  diameter 

1 

Pacific  coast 

Layer 

39 

171 

2 

Pacific  coast 

Cumulus 

44 

227 

3 

Plateau 

Layer 

30 

119 

4 

Plateau 

Cumulus 

25 

91 

5 

Eastern  U.S. 

Layer 

110 

4o4 

6 

Eastern  U.S. 

Cumulus 

4 

26 

_  1 

4o 
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TABLE  II.-  INTERVALS  OF  TEMPERATURE  AND  DROP  DIAMETER 
USED  IN  THE  DETERMINATION  OF  MARGINAL  PROBABILITIES, 
AND  MEDIAN  VALUES  IN  EACH  INTERVAL  TO  BE  USED  IN 
CALCULATING  THE  MARGINAL  LIQUID -WATER  CONTENT 


Intervals  of  temperature,  AT 
(°F) 

Intervals  of  drop  diameter,  AD 
(microns) 

Lower 

Upper 

Median  of 

Lower 

Upper 

Median  of 

limit 

limit 

. .  -  .  . 

interval,  T-j_ 

limit 

limit 

interval,  Dj_ 

20.5 

32.0 

24.0 

0 

9.5 

8 

10.5 

20.5 

15.5 

9.5 

12.5 

11 

.5 

10.5 

6.0 

12.5 

15.5 

14 

-9.5 

.5 

-4.0 

15.5 

19.5 

17 

-19.5 

1 

VO 

• 

v_n 

-13.0 

19.5 

29.5 

24 

-1*0.0 

-19.5 

-25.0 

29.5 

49.5 

37 

49.5 

00 

60 
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TABLE  III.  -  COMPARISON  OF  VALUES  OF  LIQUID -WATER  CONTENT 
FOR  PACIFIC  COAST  CUMULUS  CLOUDS,  Pe  =  O.OOl,  WITS 
CORRESPONDING  VALUES  FOR  INTERMITTENT  MAXIMUM 
CONDITIONS  LISTED  IN  NACA  TN  1855 


Temperature,  T 
(°F) 


Drop  diameter,  D 
(microns) 

Liquid  water  content,  W 
(g/m3) 

TN  1855 

Probability  analysis 

20 

2.5 

2  .46 

30 

1.3 

1.50 

50 

.4 

.^2 

20 

2.2 

2.30 

30 

1.0 

1.35 

50 

•  3 

.37 

20 

1.7 

1.45 

30 

.8 

.70 

50 

.2 

.15 

20 

1.0 

.37 

30 

.5 

.15 

50 

.1 

0 

20 

.2 

0 

30 

.1 

0 

50 

<.l 

0 
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TABLE  IV.-  COMPARISON  OF  WEIGHTED-AVERAGE  VALUES  OF 
LIQUID  WATER  CONTENT  FOR  LAYER  CLOUDS,  Pe  =  O.OOl, 
WITH  CORRESPONDING  VALUES  FOR  CONTINUOUS  MAXIMUM 
CONDITIONS  LISTED  IN  NACA  TN  1855 


Temperature,  T 
(°F) 

Drop  diameter,  D 

Liquid— water  content,  W 
(g/m3) 

(microns) 

TN  1855 

Probability  analysis 

32 

15 

0.80 

0.77 

32 

25 

.50 

.49 

32 

4o 

.15 

.20 

l4 

15 

.60 

.50 

l4 

25 

•  30 

.28 

l4 

4o 

.10 

CO 

0 

• 

-4 

15 

.30 

.17 

-4 

25 

.15 

.09 

-4 

4o 

.06 

.02 

-22 

15 

.20 

.02 

-22 

25 

.10 

.01 

-22 

4o 

.04 

.003 

-4o 

15 

.05 

0 

-4o 

25 

.03 

0 

-4o 

4o 

.01 

0 

gjggZ£S9g|j 
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TABLE  V.-  CORRESPONDENCE  BETWEEN  THE  EQUIPROBABILITY 
ANALYSIS  AND  THE  ANALYSIS  IN  NACA  TN  1855 


Type  of  cloud 
formation 

Exceed¬ 

ance 

proba¬ 

bility, 

Pe 

Approximate  corre¬ 
sponding  class  icing 
condition  as  defined 
by  TN  1855 

Equipment  to 
which  condition 
is  applicable 

Tropical  cumulus 

Pacific  coast 
cumulus 

Do  -  -  -  - 
Do  -  -  -  - 

a0.0001 

.001 

.01 

.1 

Instantaneous  maximum 
(I-M) 

None 

None 

Instantaneous  normal 
(I-N) 

Any  part  of 
airplane,  such  as 
guide  vanes,  inlet 
ducts,  where  a 
sudden  large  mass 
of  supercooled 
water  would  be 
critical,  even  for 
short  duration 

Pacific  coast 
cumulus 

Do  -  -  -  - 

Do  -  -  -  - 

.001 

.01 

.1 

Intermittent  maximum 
(II-M) 

None 

Intermittent  normal 
(II -N) 

Any  critical  com¬ 
ponent  of  the 
airplane,  where  ice 
accretions,  even 
though  slight  and 
of  short  duration, 
could  not  be 
tolerated 

Layer  clouds 
oc curring  anywhere 
in  the  U.S. (Pacific 
coast,  plateau  and 
east  coast  regions) 

.001 

.01 

.1 

Continuous  maxim-urn 
(III-M) 

None 

Continuous  normal 
(III-N) 

All  components  of 
the  airplane;  that 
is,  every  part  of 
the  airplane  should 
be  examined  with 
the  question  in 
mind,  "Will'  this  part 
be  affected  seri¬ 
ously  by  accretions 
during  continuous 
flight  in  icing 
conditions  ?" 

aExtrapolated  from  Pacific  coast  cumulus  data. 


TABLE  VI.-  EXAMPLE  PARTIAL  PROBABILITY  CALCULATIONS  FOR  PACIFIC  COAST  CUMULUS  CLOUDS 
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TABLE  VII.-  PARAMETERS  DEFINING  THE  FUNCTION,  R* 
FOR  PACIFIC  COAST  CUMULUS  CLOUDS  (CASE  2) 
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TABLE  VIII.-  EXAMPLE  CALCULATIONS  OF  EQUIPROBABILITY 
CONTOURS  FOR  SMALL  VALUES  OF  DROP  DIAMETER  FOR 
PACIFIC  COAST  CUMULUS  CLOUDS  (CASE  2) 

(a)  Pe  =  O.Ol;  T*  =  0;  Q  =  1.00 


Drop  diameter,  D 

Probability  values 

Liquid-water  content,  W 

(microns) 

S 

R 

(g/m3) 

0 

1.00 

0.0100 

1.96 

10 

.96 

.0104 

1.94 

15 

.78 

.0128 

1.88 

20 

.415 

.024 

1.67 

25 

.116 

.086 

1.24 

30 

.015 

.67 

.40 

0>)  Pe 

=  0.01;  T*  ■  32;  Q 

=  0.176 

Drop  diameter,  D 

Probability  values 

Liquid-water  content,  W 

(microns) 

S 

R 

(g/m3) 

0 

1.00 

0.057 

1.22 

10 

.96 

.059 

1.20 

15 

.78 

.073 

1.14 

20 

.415 

.137 

•  91 

25 

.116 

.490 

.4o 
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TABLE  IX.-  EXAMPLE  CALCULATIONS  OF  EQUIPROBABILITY 
CONTOURS  FOR  LARGE  VALUES  OF  DROP  DIAMETER  FOR 
PACIFIC  COAST  CUMULUS  CLOUDS  (CASE  2) 


(a)  Pe  =  O.Ol;  T»  =  0;  Q  =  1.00 


(b)  Pe  =  o.oi;  T*  =  32;  Q  =  O.176 
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TABLE  X.-  VALUES  OF  THE  FACTORS  USED  IN  THE  DETERMINATION 
OF  PARTIAL  PROBABILITIES  APPLICABLE  TO  EASTERN 
LAYER  CLOUDS  (CASE  5) 

(a)  Temperature 


Temperature 

Upper  limit  of 
tempera ture  interval 

Lower  limit  of 
temperature  interval 

Probability, 
PAT,  that 

interval,  AT 
(°F) 

Value 
of  T* 
(F°) 

Probability,  Q 
of  exceeding 

Tt 

Value 
of  T* 
(F°) 

Probability,  Q 
of  exceeding 

T* 

temperature 
lies  within 
the  interval 
AT 

32  to  20.5 

0 

1.000 

11.5 

0.520 

0.480 

20.5  to  10.5 

H.5 

.520 

21.5 

.175 

.345 

10.5  to  0.5 

21.5 

.175 

31.5 

.063 

.112 

0.5  to  -9.5 

31.5 

.063 

41.5 

.016 

.047 

-9.5  to  -19.5 

41.5 

.016 

51.5 

.003 

.013 

-19.5  to  -40.0 

51.5 

.003 

72.0 

0 

.003 

(b)  Drop  diameter 


Drop- 

diameter 

interval, 

AD, 

(microns) 


0  to  9.5  0 

9.5  to  12.5  9.5 

12.5  to  15.5  12.5 

15.5  to  19.5  15.5 

19.5  to  29.5  19.5 

29.5  to  J+9. 5  29.5 

49.5  to  00  49.5 


Lower  limit  of 


S”,  of 
exceeding 


1.000 

.780 

.480 

.230 

.090 

.022 

.0006 


Upper  limit  of 
diameter  interval 

Value 
of  D, 
(micron^ 

Probability, 
S”,  of 
exceeding  D 

9.5 

0.780 

12.5 

.480 

15.5 

.230 

19.5 

.090 

29.5 

.022 

49.5 

.0006 

00 

0 

lies  witbin  the 
interval  AD 


0.220 

.300 

.250 

.l4o 

.068 

.0214 

.0006 


H 
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TABLE  X.-  CONCLUDED 


(c)  Liquid -water  content 


Temperature 
interval,  AT 
(°F) 

Liquid -water - 
content  distribu¬ 
tion  parameters 

Liquid-water 
content,  Wp 
(g/m3) 

PWi(AT) 

^0  .63 

^0  .05 

32  to  20.5 

0.23 

0.64 

0.23 

0.63 

■  Do - 

-do- 

-do- 

.50 

•  13 

Do - 

-do- 

-do- 

.75 

.024 

Do - 

-do- 

-do- 

1.00 

.004 

20.5  to  10.5 

.24 

.72 

.24 

•  63 

Do - 

-do- 

-do- 

.50 

.18 

Do - 

-do- 

-do- 

.75 

.04 

Do - 

-do- 

-do- 

1.00 

.009 

10.5  to  0.5 

.14 

•  37 

.14 

.63 

Do - 

-do- 

-do- 

.25 

.20 

Do - 

-do- 

-do- 

.37 

.05 

Do - 

-do- 

-do- 

.50 

.009 

0.5  to  -9.5 

.10 

-=i- 

CVJ 

• 

.10 

.63 

Do - 

-do- 

-do- 

.20 

.11 

Do - 

-do- 

-do- 

.30 

.014 

Do - 

-do- 

-do- 

.4o 

.002 

-9.5  to  -19.5 

.07 

.17 

.07 

.63 

Do - 

-do- 

-do- 

.20 

.028 

Do  -  -  - 

-do- 

-do- 

•  30 

.0017 

-19.5  to  -4o.o 

.04 

0 

H 

• 

.04 

.63 

Do - 

-do- 

-do- 

.10 

.05 

Do - 

-do- 

-do- 

.15 

.004 

nggsKm 
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TABLE  XI . -  PARTIAL  PROBABILITY  VALUES  FOR  EASTERN  LAYER 
CLOUDS  (CASE  5)  FROM  DATA  PRESENTED  IN  TABLE  X 


Temperature 
interval,  AT 
(°F) 

Diameter 
interval,  AD 
(microns) 

Liquid -water 
content,  W±f 
for  AT  and  AD 
(g/m3) 

PAT 

PAD 

FWi(AT) 

APi 

32.0  to  20.5 

0  to  9*5 

0.23 

.50 

.75 

1.00 

1.30 

0.48 

0.22 

0.63 

.13 

.024 

.004 

.0005 

0.0665 

.0137 

.0025 

.0004 

.00005 

20.5  to  10.5 

0  to  9-5 

.24 

.50 

.75 

1.00 

.345 

.22 

.63 

.18 

.04 

.009 

.0478 

.0136 

.0030 

.0007 

32.0  to  20.5 

9.5  to  12.5 

.23 

.50 

.75 

1.00 

.48 

.30 

.63 

.13 

.024 

.oo4 

.0906 

.0187 

.0034 

.0006 

20.5  to  10.5 

9.5  to  12.5 

.24 
•  50 
.75 

1.00 

.345 

•30 

.63 

.18 

.04 

.009 

.0653 

.0186 

.0041 

.0009 
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Figure  2.— Graphical  presentation  of  the  equiprobability  surface  which 
presents  the  locus  of  all  combinations  of  liquid -water  content,  drop 
diameter  and  temperature  depression  below  freezing  having  the 
same  probability,  Pe,  of  being  exceeded  in  any  single  icing 

encounter. 


Liquid- water  content,  W,  g/m* 
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Figure  3 .  Constant -temperature  contours  defining  the 

equipr obabl iity  surface  for  Case  i.  (Layer  clouds 
Pacific  coast).  s 


Liquid-water  content,  W,  g/m* 


RACA  TN  2738 


-32*  F  (T'=0) 
-20*  F  (T'=  / 2) 

-  10*  F  (  T>=  22) 

-  0*F  (T'=32) 

-  iO*F  (T‘=42) 
•  20* F  (T'=52) 


20  30  40  5 

Drop  diameter,  D,  microns 

(b)  P9  *  0.0! 

Figure  J.  -  Continued. 
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Figure  3  -  Concluded , 
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NACA  TR  2738 


Drop  diameter ,  D,  microns 


(a)  Pe  =  O.i 


Figure  4.  —  Constant —temperature  contours  defining  the  equi- 
probability  surface  for  Case  2.  (Cumulus  clouds ,  Pacific 

coast). 


Liquid -water  content,  W,  g/m3 


8H 
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Drop  diameter,  D,  microns 
(b)  Pe  *  0.01 


Figure  4.  -  Continued. 


Liquid -water  content \  Wf  g/m 


58 


NACA 


Figure  4.  -  Concluded. 


Liquid  -  water  content,  Wt  g/m* 
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Drop  diameter ,  D ,  microns 


(a)  P9  r  OJ 

Figure  .5.--  Constant- temperature  contours  definina  th 
equiprobabihty  surface  for  Case  3.  (Layer  clouds 
Plateau  area).  9 


Liquid -water  content,  W,  g/m* 
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—32° 
- 20 0 
-10  0 
-  0° 
-10° 
—20° 


F  (T  s  0) 
F  (T  -  12) 
F  (T'=  22) 
F  (T\*  32) 
F  (T  -  42) 
F  (T  =  52) 


16  24  32  40 

Drop  diameter,  D ,  microns 

(C)  Pg  *  0.00/ 


Figure  5 ■  -  Concluded. 


Liquid  -water  content ,  W,  g/ ms 


62 


MCA  TET  2738 


o  10  20  30  40 

Drop  diameter ,  D ,  microns 


(a)  Pe  =  01 

Figure  6.  -  Const  ant -temperature  contours  defining  the 
equiprob ability  surface  for  Case  4.  (Cumulus 
clouds,  Plateau  area). 


Liquid -water  content,  W,  g/m s 
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Drop  diameter,  D,  microns 


(b)  Pe  r  o.O! 


Figure  6 .  —  Continued . 


Liquid-water  content W,  g/m3 
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Drop  diameter,  D,  microns 
(a)  P9  --  OJ 

Figure  7.  —  Constant— temperature  contours  defining  the 
e  quiprobabiUty  surface  for  Case  5.  ,  (Layer  clouds. 
Eastern  United  States). 


Liquid- water  content ,  Wt  g/m 
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Drop  diameter,  D,  microns 


(b)  Pe  --  0.0/ 


Figure  7.  -  Continued. 


Liquid -water  content,  W,  g/m 
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8  16  24  32  40  48 

Drop  diameter,  D,  microns 

(c)  p9  -  0.00/  -  y  - 


Figure  7.  -  Concluded. 
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Figure  8-  Variation,  with  horizontal  extent,  of  the  factor  by  which  values  of  liquid-  water 
content  from  the  equiprobability  charts  should  be  multiplied  to  include  the  effect 
of  horizontal  cloud  extent  for  different  exceedance  probabilities. 
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Liquid -water  content  W/t  grams  per  cubic  meter 


(a)  Case  /  -  Layer  clouds ,  Pacific  Coast 
Figure  10.-  Partial  probability  charts. 


Partial  probability  APj,  dimensionless 
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Liquid -water  content  Wf,  grams  per  cubic  meter 
0  .5  i.O  1.5  2.0  2.5  3.0  3.5 


(b)  Case  2  -  Cumulus  clouds,  Pacific  Coast. 
Figure  10.  — Continued. 
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fo/  Case  3  -  Layer  clouds ,  Plateau  Region. 


Figure  !0.  —  Continued. 
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Liquid -water  content  Wj,  grams  per  cubic  meter 


(d)  Case  4  -  Cumulus  clouds,  Plateau  Region. 
Figure  10.-  Continued. 


Partial  probability  APj,  dimensionless 
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Liquid -water  content  Wj,  grams  per  cubic  meter 


(e)  Case  5  —  Layer  clouds ,  Eastern  United  States 
Figure  10.  —  Concluded. 
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Probability,  R,  dimensionless 


S9  .90  70  .50  .30  JO  .05  .02  .0/  .005 


Percent  of  icing  encounters  with  maximum  liquid  water  less  than  W 


Figure  12 -Frequency  distributions  of  maximum  liquid -water  content  per 
icing  encounter  for  two  values  of  maximum  temperature  for 
Pacific  Coast  cumulus  clouds. 
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Figure  13-  Distribution  parameters  which  define  the  probability 
function  R(W,  T'J  for  Pacific  Coast  cumulus  clouds. 
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Cumulative  frequency,  percent  "\naca^ 

Figure  15.  -  Frequency  distribution  of  maximum  value  of  mean- effective 
drop  diameter  observed  in  each  icing  encounter  for  Pacific  Coast 
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Drop  diameter,  D,  microns 


Figure  17. — Parameters  defining  the  distribution  of  liquid 
water  content  for  values  of  drop  diameter  greater  than 
D  for  Pacific  Coast  cumulus  clouds ;  T‘  A  0. 
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Drop  diameter ,  D,  microns 


Figure  18.-  Equiprobability  contours  constructed 
from  the  data  of  tables  VUL  and  IF ; 


Pe  =  O.Ot. 
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freezing  of  the  most  severe  icing  condition  per  icing 
encounter  for  Eastern  U.  S.  layer  clouds. 


Liquid  water  content,  W,  9/m3 
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Probability,  R,  dimensionless 


Figure  20.— Frequency  distribution  of  the  maximum  liquid -water 
content  observed  during  each  icing  encounter  with  T'  equal 
to  or  greater  than  the  indicated  value. 
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Probability,  S,  dimensionless 
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Probability ,  5/  dimensionless 


Percent  of  icing  encounters  with  maximum  observed  value  of  drop 

diameter  less  than  D 


Figure  23.-  Frequency  distribution  of  the  largest  value  of  drop  diameter 
observed  during  each  icing  encounter ;  layer  clouds.  Eastern  United 

States. 
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Drop  diameter ;  D,  microns 


Figure  24.  —  Distribution  parameters  defining  the  function  R'(W,D) 
for  the  case  of  T'  0;  layer  clouds ,  Eastern 

United  States. 
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Probability,  S,“  dimension/ess 

10  .05  .02  .01  .005  .002  .001  .0005 


Figure  26- Composite  frequency  distribution  of  drop  diameter  used  in  calculating  partial  probabilities  for 

layer  clouds  in  Eastern  United  States. 
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Log  PAT  +  Log  PWj(AT) 


Figure  27  —  Example  construction  of  a  partial 

probability  chart 


Liquid-  water  content ,  IV j,  grams  per  cubic  meter 
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Liquid -water  content,  W/,  grams  per  cubic  meter 


(b)  Step  2 


Figure  27.  —  Continued. 
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Liquid-  wafer  content ,  W/,  grams  per  cubic  meter 
0  25  .50  75  LOO 


Q> 

£ 

>4 


(c)  Step  3 

Figure  27  -  Concluded. 
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